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ABSTRACT 

C*") \ We use hydrodynamical models of the wind-collision region (WCR) in the archetype 

colliding-wind system WR 140 to determine the spatial and spectral distribution of the 
' radio, X-ray and 7-ray emission from shock accelerated electrons. Our calculations are 

for orbital phase 0.837 when the observed radio emission is close to maximum. Us- 
ing the observed thermal X-ray emission at this phase in conjunction with the radio 
emission to constrain the mass-loss rates, we find that the O-star mass-loss rate is con- 
sistent with the reduced estimates for 04-5 supcrgiants bv lFullerton. Massa fc Prinial 
^5 ' l)2006j) . and the wind momentum ratio, 77 = 0.02. This is independent of the opening 

\Q \ angle deduced from radio VLBI observations of the WCR that we demonstrate fail to 

C^) ■ constrain the opening angle. 

We show that the turnover at ~ 3 GHz in the radio emission is due to free- 
free absorption, since models based on the Razin effect have an unacceptably large 
fraction of energy in non-thermal electrons. We find the spectral index of the non- 
O [ thermal electron energy distribution is flatter than the canonical value for diffusive 

shock acceleration (DSA), namely p < 2. Several mechanisms are discussed that could 
' lead to such an index. Our inability to obtain fits to the radio data with p > 2 does 

not exclude the possibility of shock modification, but stronger evidence than currently 
exists is necessary for its support. 

Tighter constraints on p and the nature of the shocks in WR 140 will be obtained 
from future observations at MeV and GeV energies, for which we generally predict 
■ lower fluxes than previous work. Since the high stellar photon fluxes prevent the 

acceleration of electrons beyond 7 > 10 5 — 10 6 , TeV emission from colliding-wind 
binary (CWB) systems will provide unambiguous evidence of pion-decay emission 
from accelerated ions. We finish by commenting on the emission and physics of the 
multiple wind collisions in dense stellar clusters, paying particular attention to the 
Galactic Center. 

Key words: stars:binaries:gcncral - stars:carly-typc stars:individual:WR 140 - 
stars:Wolf-Rayet - radio continuumrstars - Xraysrstars 
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1 INTRODUCTION 

Massive binary systems containing a Wolf-Rayet (WR) 
and O-type star often exhibit synchrotron emission aris- 
ing fr om relativistic electrons in th e presence of a magnetic 
field jDoughertv fc Williams! l200fA . The non-thermal elec- 
trons are widely thought to be accelerated through diffu- 
sive shock acclerati on (DSA), though other mechanisms are 
also possible (e.g., IJokipiil |l987l : Ijardine. Allen fc Pollock! 
Il996f) . High-spatial-resolution radio observations (e.g., 
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Williams ct al. lll997l:lpoughertv. Williams fc Pollaccoll2000l: 
Dougherty et alj|2005l) " have revealed that the acceleration 
site is the region where the massive stellar winds collide - the 
wind-collision region (WCR). Thus, colliding-wind binary 
(CWB) systems present an important laboratory for inves- 
tigating the underlying physics of particle acceleration since 
they provide access to higher mass, radiation and magnetic 
field energy densities than in supernova remnants (SNRs), 
which have been widely used for such work. 

CWB systems are also ideal for understanding related, 
but more complicated, systems involving colliding winds, 
such as the cluster of massive stars in the central par- 
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sec of the Galaxy jGhez et alJ 12005ft . and young, dense 
stellar clusters (s uper star clust ers), such as the A rches 
llFieer et alJ Il999l). NGC 3603 llMoffat et alJ 12(1 



.... , UMmn and 

R136 jFiger. McLean fc Morrisl 119931 . Non-thermal radio 
emi ssion has recently been detected from the Arches clus- 
ter ijYusef-Zadeh et alll2003fl , and the non-thermal radio fil- 
aments and "streaks" near the Galacti c Center may hav e 
their origin in CWBs and cluster winds jYusef-Zadehl2003l) . 

In this paper we model the thermal and non-thermal 
radio, X-ray and 7-ray emission from one of the most exten- 
sively observed CWB systems, WR140. In Sec.[5]we describe 
our assumptions for the non-thermal particle spectrum. 
Sec. |3 contains details of our implementation of the various 
processes responsible for the non-thermal X-ray and 7-ray 
emission. Details of the modelling of the radio emission, in- 
cluding several cooling processes, and vari ous emission and 
absorp tion mechanisms, c a n be found in iDoughertv et all 
J2003T) and iPittard et all fcOOrl) . The application of our 
model to observations of WR140 is described in Sec. [I] and 
the implications of our modelling are discussed in Sec.|^| In 
Scc.|S] wc summarize and note future directions. 



2 THE NON-THERMAL PARTICLE 

SPECTRUM AND MAGNETIC FIELD 

In this work, it is assumed that the non-thermal particle dis- 
tribution can be approximated as a power-law n(-y) oc 7~ p , 
where 7 is the Lorentz factor and p is a free-parameter de- 
termined from model fits (primarily) to the radio emission. 
A number of mechanisms can produce such a distribution, 
as we now discuss along with our assumptions concerning 
the normalization of the non-thermal particle spectrum and 
the magnetic field within the WCR. 

2.1 Acceleration mechanisms and spectral index 
of non-thermal particles 

There are several possible mechanisms for accelerating par- 
ticles in CWB systems, including DSA, reconnection, and 
many turbulent processes. The spectral index of the non- 
thermal particles can, in each instance, span a wide range of 
possibilities, as discussed below. An estimate of which mech- 
anism is likely to be dominant can be made by comparing 
their relative efficiencies. 



2.1.1 Reconnection 

In CWB systems, particles may be accelerated at the con- 
tact discontinuity between the two winds as magnetic field 
lines from the two stars are for ced together and reconnect 
jjardine. Allen fc Pollocklll996l) . and/or within the volume 
of the WCR as tangled field lines reconnect. Reconnection 
can occur either as a quasi-steady, or as an explosive, pro- 
cess. The energy released is transferred into thermal energy 
through direct Joule heating of the plasma, into the kinetic 
energy of macroscopic motions, and into accelerated parti- 
cles. Unfortunately, the role of magnetic reconnection in ac- 
celerating particles is complex, not least because of the many 
different types of reconnection, and the intricate magnetic 
environment that is usually associated with reconnection. 
This environment can accelerate particles directly (through 



strong electric fields), stochastically (i.e. through the second- 
order Fermi process due to high turbulence), and at the 
MHD shock waves which are an int egral part of a recon- 
necting field (i.e. through DSA - see iPriest fc Forbesllioobl 
and references therein). Reconnection is usually thought to 
produce steep spectra i.e. p > 2, such as is seen in solar 
flares. 

A rough upper limit of the power available for particle 
acceleration due to reconnection can be estimated as follows. 
If reconnection occurs as a quasi-steady process, the power 
released in a current sheet of area A is 
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where R m is the magnetic Reynold's number, B is the mag- 
netic field strength, and va is the Alfven speed llPetschekl 
Il964l) . If reconnection occurs at the contact discontinuity, 
the area of the current sheet, A, is approximately the area of 
the apex of the WCR i.e. A ~ -nr Q , where r = r/ 1/2 D/(l + 
77 1 / 2 ) is the distance from the stagnation point of the WCR 
to the center of the O star, D is the separation of the stars, 
and the wind momentum ratio 77 = Movo/-MwrVwr where 
Mo and M^h are the mass-loss rates of the O and WR stars 
with terminal wind speeds of vo and «wr respectively. For 
our model B (see Table [TJ, R m ~ 10 20 , A ~ 5 x 10 27 cm 2 , 
v A ~ 6xl0 7 cm s" 1 , and B = 1.2 G, so P ~ 3xl0 32 ergs" 1 . 
Since P is strongly dependent on B, the available power from 
reconnection could be considerably higher than this esti- 
mate, and may be further augmented by reconnection occur- 
ing throughout the WCR region. While we expect that only 
a fraction of this total power will end up in non-thermal par- 
ticles, there is probably sufficient reconnection luminosity 
to produce the synchrotron emission observed from WR 140 
(~ 5 x 10 32 ergs -1 ). However, it may struggle to supply the 
energy in IC 7-ray emission, which is expected to be ~ 10 x 
greater (see Table 0. 



2.1.2 Turbulent mechanisms 

Particles may also be accelerated stochastical ly, by magnetic 
scattering off turbulent motions ([Fermi 1949). In this mech- 
anism, non-thermal particles either gain or lose energy de- 
pending on whether the magnetic field is moving towards 
or away from the particle. However, as the scattering prob- 
ability is proportional to the relative velocity between the 
particle and the plasma, there is a slight excess of head-on 
scatterings and a mean increase in momentum. This pro- 
cess is known as the second-order Fermi mechanism, and in 
principle, hard power-law spectra (p < 2) may be obtained. 
The second-order Fermi mechanism is slower and less effi- 
cient than DSA (a first-order Fermi process) in our mod- 
els of WR 140, since the pre-shock Alfvenic Mach number, 
Ma = v^Anp/B > 3 - see Tabled Another possibile mech- 
anism is acceleration in turbu lence from resistiv e magne- 
tohydrodynamical fluctuations ijNodes et al.ll20f3) that can 
produce hard power-law spectra with p < 2. 



2.1.3 Diffusive shock acceleration 

DSA is a variant of the original Fermic mechanism, in which 
particles are accelerated by scattering back and forth across 
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a shock. The key difference is that fast particies gain en- 
ergy each time they are scattered back across the shock, 
until their eventual loss downstream of the shock. To a first 
approximation, DSA produces a p ower-law energy dis tribu- 
tion of non-thermal particles fsee lPittard et afll2006t for a 
discussion of departures from this). The maximum power 
that is available to DSA is the kinetic luminosity of that 
part of the stellar winds which are processed through the 
WCR (in model B this is ~ 10 36 ergs -1 ). Therefore, DSA 
easily meets the energy requirements for producing the syn- 
chrotron emission seen from WR140 (~ 5 x 10 32 ergs -1 ). 

At a single, parallel shock the value of p is determined 
by the compre ssion ratio of the scattering centers across the 
shock, rk (e.g., Bell 1978). In the strong shock limit, the gas 
compression ratio r — 4; if « r, then p ~ 2. However, 
p can be greater or smaller than this value depending on 
the values of r and the ratio of gas to magnetic pressure in 
the u pstream flow, which together determine iSchlickeiserl 
120021) . 

There are additional ways in which the spectral in- 
dex can vary from the canonical strong shock limit of 
p = 2. For instance, when particle acceleration is efficient, 
the diffusion of non-thermal particles upstream of the sub- 
shock exerts a significant back-pressure on the pre-shock 
flow, causing the gas velocity to decrease prior to the sub- 
shock in a smooth precursor (Fig. Q. As the highest en- 
ergy particles are able to diffuse further upstream from 
the subshock, and see the greatest velocity difference be- 
tween the upstream and downstream flow, a concave cur- 
vature to the non-thermal particle energy spectrum is cre- 
ated. In addition, the velocity jump across the subshock 
in this situation is much less than in the classical shock, 
which results in a reduced post-shock thermal tempera- 
ture. Convincing evidence for such "shock modification" 
is provided by analysis of observations o f several young 
SNRs, including Tycho JVolk et alj l2002Tl and SN 100 6 
feamba et alj 120031: iBerezhko. Ksenofontov fc Volkl Exit) . 
and from direct spacecraft measurements of the Earth's 
bowshock and interplanetary shocks (see references in 
Ieu ison fc Cassam-Chenaill2005f) . Non-thermal ions are be- 
lieved to be primarily responsible for modifying shocks be- 
cause they are thought to dominate the non-thermal energy 
density (see Sec. 12. 21 . Since the diffusion length of electrons 
is much smaller than ions for a given Lorentz factor, non- 
thermal electrons are confined to a region much closer to the 
subshock. Thus, only those electrons with momenta > m p c 
(i.e. with Lorentz factors > 2000) are able to diffuse far 
enough upstream to see the velocity difference required to 
obtain a spectral index p < 2 (assuming p = 2 for the un- 
modi fied shock - see Fig. 1 in lEllison. Decourchelle fc Balled 
120041) . In CWB systems, the synchrotron emission at GHz 
frequencies is from electrons with relatively low momentum. 
For example, at the apex of the WCR in model B of this 
work, electrons with 7 = 50 (i.e. momentum p = 50 m c c ~ 
0.03 m p c) produce synchrotron emission which peaks at 
v = 0.29^ c ~ 3.6 GHz, where v c is the frequency where the 
synchrotron spectrum cuts off iPacholczv kjll97pj). Hence, if 
shock modification occurs in these systems, the synchrotron 
emission should be characterized by p > 2, though if re- 
acceleration occurs this is not necessarily the case. 

The spectral index of non-thermal particles can be mod- 
ified if they pass through a sequence of shocks. In the limit 



of many shocks, the spectral index approaches p = 1, in- 
dependent of the sho ck strength. This proc ess is known as 
"re-acceleration" fsee lPope fc MelroseHl994l . and references 
therein), and if the particles suffer energy losses between the 
accelerating shocks, p can bec ome even smaller fechlickeiserl 
ll984HMelrose fc Croucbll997D . Shock modification and par- 
ticle re-acceleration are not necess arily mutually ex cl usive. 

Some previous works llKirk fc Heavensl Il989l : 
iNaito fc Takaharal | l995f ) have also claimed that p < 2 
can be obtained if the shock is oblique (i.e. when the B-field 
is at some angle to the shock normal), but these rely on 
a non-standard shock acceleration model, such as non- 
isotropic scattering, and we do not consider this further. 
Mechanisms which increase the effect ive Mach number o f 
the shock may also yield p < 2 (e.g., iMalkov et al] 1200(1) . 
while p > 2 may be obtained if sub-diffusive cross-fiel d 
transport is important jPuffv et alJll995l : iKirk et alJll996tl . 

2.1.4 Summary of assumed NT particle spectrum 

From energy considerations, DSA and reconnection appear 
to be the most efficien t mechanisms, w ith DSA likely hav- 
ing the edge (see also lBednareklE005f) . Therefore, for the 
remainder of this work we assume that the non-thermal 
particles (electrons and ions) are created by DSA at the 
global shocks bounding the WCR and have power-law en- 
ergy distributions. We further assume that the downstream 
magnetic field is highly tangled. While we do not model the 
actual acceleration process itself, this is consistent with the 
Bohm limit where the magnetic fluctuations are of the same 
strength as the background field i.e. SB/B ~ 1. The strong 
cross-field diffusion which results in this limit allows efficient 
injection at quasi-perpendicular shocks. 

The diffusion length scale of accelerating particles is 
given by idtff = ft/ushk, where u s hk is the shock speed, 
K — Xv/3 is the spatial diffusion coefficient, A ~ 3r g is 
the mean-free-path, and r g = 7771c 2 jqB is the gyroradius 
of particles with Lorentz factor 7, mass m, and charge q in 
a magnetic field B. The acceleration of non-thermal elec- 
trons at the shocks is not affected by the flow geometry of 
the downstream plasma, as the diffusion lengthscale is only a 
small fraction of the width of the WCR (e.g., l^is ~ 10 10 cm 
for electrons with 7 = 10 5 at the shock apex of model B, 
whereas the distance between the global shocks and the con- 
tact discontinuity is ~ 2x 10 13 cm). Therefore, the geometry 
of the post-shock flow will have no impact on our calculation 
of the synchrotron and IC emission, or our finding that the 
electron spectrum has p < 2. The flow geometry may be- 
come important for non-thermal ions, but for protons and 
He nuclei this is only likely to be the case for 7 > 10 4 . De- 
tailed work outside the scope of this paper is required to 
examine the effect of the flow geometry on our predictions 
of the high energy emission from the decay of neutral pions. 

We do not consider the inverse-Compton (IC) cascade 
of secondary electrons and positrons that occur when the 
optical depth of 7-rays to pair-production in a stellar radia- 
tion field is high and when IC scattering is efficient, though 
in WR 140 the secondary 7-rays produced in the cascade and 
which subsequently escape the system may be comparable 
to the observed flux of primary 7-rays above Ge V energies 
for th e orbital phase considered in this work (cf. iBednarekl 
120061) . The non-thermal electrons formed by hadronic colli- 
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Figure 1. Top: Schematic structure of a strong hydrodynamic 
shock. If DSA places a negligible amount of energy into acceler- 
ated particles, an ordinary discontinuous classical shock occurs 
(dashed line). In this case the undisturbed prc-shock flow is sep- 
arated from the post-shock flow by a discontinuity known as a 
subshock across which the flow velocity decreases, and the den- 
sity and pressure increase. The scale length of the subshock, Z s , 
is determined by microphysical dissipation processes that result 
in gas heating. In contrast, when particle acceleration is efficient, 
the diffusion of non-thermal particles upstream of the subshock 
exerts a back-pressure on the pre-shock flow, causing the gas 
velocity to decrease prior to the subshock in a smooth precur- 
sor (solid line). The scale length of the precursor, (>> l B ), 
is approximately the diffusion length of the highest energy ions. 
Middle: The spatial distribution of the number density of non- 
thermal particles as a function of energy in the modified shock. 
The highest energy particles are able to diffuse further upstream 
from the subshock, and see the greatest velocity difference be- 
tween the upstream and downstream flow. Non-thermal particles 
in the precursor are typically scattered back across the subshock, 
and either undergo further acceleration, or escape downstream. 
Bottom: The non-thermal particle energy spectrum multiplied by 
7 2 to emphasize the spectral curvature of the non-linear (i.e. effi- 
cient particle acceleration) case. When DSA places relatively little 
energy into accelerated particles the standard test particle spec- 
trum is obtained, as shown here for p = 2 (solid). In contrast, the 
modified shock produced when DSA is efficient leads to a concave 
cur vature in the non-thermal particle en ergy spectrum (dashed - 
see Ell ison. Dccourchcllc & Ballet 20o3i and references therein). 



sions, and the knock-on electrons formed when non-thermal 
nuclei collide with thermal (atomic and free) electrons, are 
also ignored in this work, though the emission resulting from 
these processes is not expected to be dominant (see also 
| MastichiadiJl996tlBaring et alJl999HBerrington fc Dermerl 
120031) . We also ignore an additional (though likely small) 
contribution to the non-thermal radio emission from parti- 
cles accelerated at wind-embedded shocks caused by the line- 
deshadowing instability in radiatively driven winds, since in 
WR 140 the non-thermal emission is clearly related to bina- 
rity. We set p as a free parameter for the particle distribution 
immediately downstream of the shocks, with the assumption 
that it is spatially invariant along the face of the shocks, 
and identical for electrons and ions. For the ions p does not 
change in the post-shock flow, but as the non-thermal elec- 
trons advect downstream IC cooling modifies their energy 
distribution, and in this case p is dependen t on both their 
spatial position and 7 iPittard et alj|200fj) . For instance, 
in model B, electrons accelerated at the WR shock on the 
line of centers between the stars cool to 7 < 100 within a 
time of « 6.3 x 10 4 s, during which they are advected a dis- 
tance of « 4.7 x 10 12 cm (0.015 D) downstream (assuming 
that the speed of the post-shock flow is 0.25 Voo in this re- 
gion of the WCR). Electrons accelerated at a point on the 
WR and O shocks which is a distance D from the O star 
cool to 7 < 1000 within a time of ~ 2 x 10 5 s. During this 
time they are advected a distance of 5.9 x 10 13 cm (0.19 D) 
downstream (assuming a post-shock flow speed close to the 
terminal wind speeds). 



2.2 Normalization 

2.2.1 The non-thermal particle spectrum 

Since the population of non-thermal particles is not deter- 
mined from first principles in our models, the non-thermal 
electron energy density, U te \, e , is assumed to be propor- 
tional to the thermal particle internal energy density, [/th 
i.e. {7 re i,e = Croi, c tAh- Also, it is assumed that shock obliq- 
uity does not influence the normalization of the non-thermal 
electron energy density. However, we note that in gen- 
eral, the higher the obliquity, the lower the injection ef- 
ficiency (although particles are actually accelerated more 
rapidly as a result of shock drift along the surface and 
slower diffusion in the shock normal direction - see, e.g., 
[ Ellison. Baring fc Jonesll995l.ll99t3:lEllison. Jones fc Baring! 
1999). In wide binaries, such as WR 140, the shocks near the 
head of the WCR should, in principle, be highly oblique. 
Yet this picture is further complicated if the level of mag- 
netic turbulence is at the Bohm limit, since then the shock 
does not have a well-defined obliquity, and strong cross- 
field diffusion results in effic ient injection (see Fig. 6 in 
lEllison. Baring fc Jonesl Il995h . Such concerns are beyond 
the scope of the present work. 

To obtain an estimate of the 7-ray flux from the de- 
cay of neutral pions, the energy density of relativistic ions 
must be specified. Similar to the non-thermal electrons, it 
is assumed that {7 r8 i,i = Crei.i^th- Fits to the radio data are 
only sensitive to £rei,e and constraining £ rc i.i requires unam- 
biguous detection of the emission from neutral pion-decay, 
which has not been attained for CWB systems. Moreover, 
the relative injection efficiency of electrons and ions is not 
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known, while the value of Gei.i/Cre i.e also depends on the 
spectral index p iSchlickeiserll2002l) .' Having said this, it is 
reasonably certain that DSA transfers considerably less en- 
ergy into non-thermal el ectrons than into non-thermal ions 
(e.g.. iBaring et al.lll999h . while the relative difficulty of ac- 
celerating electrons is also consistent with the well estab- 
lished value for the ratio of proton to electron energ y densi- 
ties for Galactic cosmic rays of ~ 100 llLongairll994l) . In this 
work we somewhat arbitrarily assume that Ciei,i = 100£ re i,e, 
unless otherwise noted, though the actual ratio in WR 140 
could differ significantly. As a result, our predictions of the 
7-ray flux from neutral pion decay and neutrino flux from 
charged pion decay are highly uncertain. However, hard up- 
per limits for these fluxes can be obtained due to the re- 
quirement that the total power put into non-thermal par- 
ticles must be less than the available kinetic power in the 

WCR i.e. Crel.e + Crel.i < 1. 

As non-thermal particles flow downstream, they are 
subject to adiabatic losses. This is achieved in an approxi- 
mate fashion in this work through the assumption that U Te i,e 
and E/ re i,i are proportional to U t h- In reality the adiabatic 
indices for thermal and relativistic plasma differ, and adia- 
batic expansion also reduces the Lorentz factor of the parti- 
cles. While neither of these effects is treated properly in our 
model, the specific parameterization we have chosen will not 
lead to significant errors in the resulting non-thermal radio 
emission. This is because IC cooling is so rapid for high en- 
ergy electrons that synchrotron emission above ~ 1 GHz is 
largely confined to a thin region near the shocks, where the 
effects of adiabatic expansion of the flow are minimal. Like- 
wise, the radiation from pion decay is concentrated near the 
apex of the WCR. 



2.2.2 The magnetic field 

If the field is frozen into the post-shock plasma of density 
p, one would anticipate B oc p (Ub oc p 2 ), but amplific ation 
of th e magnetic field by the accelerated particles (e.g.. iBell 
l2004f) and/or reconnection complicates this simple scaling. 
This also implies that it is not necessarily straightforward to 
determine the B-field in the WCR from assumptions about 
the strength of the surface magnetic field at each star. In- 
direct evidence for a mplification has recently been obtained 
in several SNRs (e.g., Berezhko. Ksenofontov fc Vol]Jl2003t 
iBerezhko fc Volklbooi lEllison fc Cassam-Chena]||2005D . In- 
stead, for convenience, we have assumed that the magnetic 
energy density is proportional to the thermal particle energy 
density i.e. Ub = CstAh- This choice is equivalent to some 
amplification in the downstream flow, but we have verified 
that there is little difference in the shape of the synchrotron 
spectrum between these two methods. Compared to other 
work where the B-field is assumed to be spatially invari- 
ant, this approximation at least has the benefit that the B- 
field declines away from the symmetry axis. Undoubtedly, a 
better description for the B-field will be obtained in future 
MHD models of the WCR. Our final assumption is that the 
post-shock B-field is highly tangled so that the synchrotron 
emission is isotropic. 



2.2.3 Limits on Cb and Crei.e 

When modelling the radio emission from specific systems, 
Cb and Crel.e are chosen to best match the observed radio 
emission. As Cb decreases, the value of Crel.e has to increase 
to reproduce an observed radio flux. A lower limit on Cb can 
be obtained from the fact that the Razin turnover frequency, 
!/r, increases with declining Cb (^r ~ 20n c /B, where n c is 
the electron number density and the magnetic field, B = 
\/87t{/b), or because the value of Crei.e required to match 
an observed flux exceeds the available power for particle 
acceleration in the WCR. An upper limit on Cb is obtained 
from the necessary c onditions for particle acceleration (see 
lEichler fc Usovll993D . Since there is currently little guidance 
for the ratio of CB/Crei,e in CWB systems, we consider this 
to be a free parameter in our models, subject to the above 
constraints. 



2.3 Maximum energy of non-thermal particles 

The maximum value of 7 for electrons is obtained by balanc- 
ing the rate of energy increase from DSA to t he rate of en- 
ergy loss by IC cooling (e.g., see Eqs. 4 and 7 in lPittard et alJ 
2006). This is spatially dependent due to the geometric di- 
lution of both the stellar radiation and magnetic fields, and 
also dependent on orbital phase for binaries with eccentric 
orbits like WR 140. An evaluation of these equations at the 
points where the shocks intersect the line of centers through 
the stars reveals that 7 ma x ~ 10 at orbital phase 0.837 for 
sensible values of B. Most of the wind power passing into the 
WCR occurs wi thin an off-axis distance of ~ 2ro from the 
line of centers dPittard fc Stevens] l2002f) . and one expects 
that most of the energy transfer into non-thermal particles 
must occur within this volume too. The value of 7 max eval- 
uated at a position on the shock front at the extremity of 
this region is only 25 per cent higher than on the line of 
centers if the magnetic field strength is assumed to decline 
as 1/r and if there is no dependence on the shock obliq- 
uity (as already mentioned, this is not well defined if the 
level of magnetic turbulence is at the Bohm limit). There- 
fore, for simplicity, throughout this work it is assumed that 
7max (electrons) = 10 along the entire shock fronts. 

While the above evaluation was for the Thomson limit, 
at such energies IC losses occur in the Klein-Nishina regime, 
and the spectral distribution as well as the cut-off energy will 
change. However, in the calculations in this paper, IC cool- 
ing is so rapid near the apex of the WCR that the down- 
stream population of non-thermal electrons in this region 
is essentially independent of whether Thomson or Klein- 
Nishina cross-sections are used. Thus, the IC emission and 
our predictions for the 7-ray fluxes are essentially identical 
too, with the exception that the high energy cutoff is sensi- 
tive (by a factor of order 2) to such details. 

The maximum non-thermal ion energy is set by bal- 
ancing the timescale of their accelerati on to their advec tion 
timescale out of the system. Following iBednareld J2005I) 1 . 

,. . 3ZeBr v 
7max(l ° ns) = (2) 

1 In Eq. 3 of Bcdnaiej]| l2005l> the denominator should have a 
factor c 3 as shown here in Eq. [2] 
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where Ze and Am p are the charge and mass of the ion, and B 
and v are the pre- shock magnetic field and wind speed. For 
wide binaries where the pre-shock magnetic field declines as 
1/r, 7max(ions) is expected to be independent of D. 

For the models considered in this paper, 7 max (ions) ~ 
10 , so for convenience, we set 7 max = 10 5 for both the 
electron and ion distributions in the remainder of this work. 
The predicted 7-ray flux is, of course, sensitive to the value 
of 7max chosen for both the electrons and ions. In addition, 
7min = 1 is assumed. 

2.4 Composition of non-thermal nuclei 

In the simplest test-particle calculations of DSA i.e. with no 
feedback, all particles accelerated in a given shock will have 
identical power-law spectral shapes in momentum, regard- 
less of their charge, and the non-thermal particle compo- 
sition reflects the composition of thermal particles. How- 
ever, in non-linear models, the spectral shape and injec- 
tion and acceleration efficiencies are dependent on the ionic 
species. Ions with large mass-to-charge (A/Q) ratios diffuse 
further upstream than protons of the same energy per nu- 
cleon (provided that both are non-relativistic) , and are ac- 
celerated more efficiently and obtain a flatter spectrum be- 
cause they "see" a larger velocity difference and are more 
easily injected to suprathermal energies. Among the volatile 
elements, observational support comes from Galactic cos- 
mic rays, which show an enhancement of heavier elem ents 
relative to lighter ones (iMever. Drurv fc Ellison! Il997t) . A 
general enhancement of the refractory elements relative to 
the volatile ones i s also apparent from the Gala ctic cosmic 
ray composition iMever. Drurv fc Ellison! Il997l) . This has 
been explained by acceleration of dust grains (which can 
have huge mass-to-charge ratios of ~ 10 4 — 10 8 ), followed 
by the sputtering of atoms from these grains. However, the 
harsh UV radiation fields in CWB systems makes it unlikely 
that dust grains can form near the head of the WCR, so an 
enhancement of refactory elements is not expected. Several 
CWB systems, including WR 140, show observational signa- 
tures of either episodic or continuous dust formation (e.g ., 
IWilliamslll99r3 : iTuthill et al"1ll999l: iMarchenko et al.ll20021) . 
though in all cases the dust is thought to form well down- 
stream of the apex of the WCR. Observations have further 
revealed that th e relative abundance of heavy ions inc reases 
with energy (see lBvkov fc ToptygirJl200ll ; lHillasll2005l) . This 
is supported by the results from non-linear models where 
variation with shock speed an d Mach number is also seen 
jEllison. Drurv fc Meverll997l) . It is also possible that shock 
obliquity may influence the non-thermal particle composi- 
tion. 

High energy heavy nuclei may also lose nucleons due to 
collisions with stellar photons i.e. by photo-disintegration. 
This occurs if the photon energy in the reference fr ame of the 
nucle i, E 1 = 3£;bT7(1 + cos 8), exceeds ~ 2 MeV iBednarekl 
I2005T) . Here T is the temperature of the radiation field. In 
WR 140, this process applies to nuclei with 7 > 8.5 x 10 4 , and 
is only likely to affect a small part of the top end of the non- 
thermal particle distribution. In addition, an estimate of the 
ratio of the characteristic photo-disintegration timescale to 
the acceleration and advection timescales reveals that photo- 
disintegration is not efficient for particles with 7 < 10 5 in 
WR 140 at the orbital phase considered in this work. The 



fragmentation of high energy nuclei through collisions with 
thermal ions is likewise negligible. 

Given the uncertainties discussed earlier, we adopt the 
following simplifications: i) identical spectra for all non- 
thermal particles, ii) identical composition to that of the 
thermal particles, iii) no dependence on the shock obliquity. 



3 NON-THERMAL X-RAY AND 7-RAY 
EMISSION AND ABSORPTION 

The presence of non-thermal electrons and ions should give 
rise to X-ray and 7-ray emission from several separate mech- 
anisms, including IC scattering, relativistic bremsstrahlung, 
and pion decay (the non-thermal electrons are not ener- 
getic enough to produce synchrotron X-ray emission). While 
definitive evidence for non-thermal X-ray and 7-ray emission 
from CWBs does not yet exist, many of the yet unidenti- 
fied EGRET s ources appear correlated with pop ulations of 
massive stars iRomero. Benaglia fc TorresHl999l) . and some 
systems appear to show a power-law tail e xtending to higher 
energies than the therm al X-ray emission JViotti et all2004l : 
iDe Becker et alJl2004bT) . A critical question is whether the 
non-thermal X-ray and 7-ray fluxes from CWB systems are 
sufficiently high to be detected with instruments onboard 
current observatories, and instruments on observatories that 
will be operational in the near future. 

Previous predictions of the IC emission have usually 
been based on the fact that the ratio of the luminosity from 
IC scattering to the synchrotron luminosity is equal to the 
ratio of the photon energy density, f7 p h, to the magnetic field 
energy density, Ub' 

Lie _ Uph 
Lsync Ub 

Typically, the ratio of Uph/Us is evaluated only at the stag- 
nation point, whereas in reality the emission is generated 
throughout the WCR, with Uph/Us varying spatially. In ad- 
dition, L sync is usually set to the observed synchrotron lu- 
minosity, whereas free-free absorption by the extended win d 
envelopes can be significant (see, e.g.. IPittard et al"ll2006t) . 
In such cases the intrinsic synchroton luminosity, and con- 
sequently the non-thermal X-ray and 7-ray luminosity, will 
be underestimated if radio data, particularly at lower fre- 
quencies, are used without regard to free-free absorption (cf. 
Figs. 1141 and 1151 . Furthermore, the predictive power of 
Eq. [3] is greatly undermined by the fact that the magnetic 
field strength in the WCR is generally not known with any 
certainty. Since Ub oc B 2 , small changes in the estimated 
value of the magnetic field stength, B, can lead to large 
changes in L ic (see Fig. 3 in lBenaglia fc Romeroll2003t) . 

To obtain a more robust estimate of the IC emission 
from CWB systems, the population and spatial distribution 
of non-thermal electrons must be determined directly. This 
information can, in principle, be acquired from model fits 
to radio data, since the IC emission arises from the same 
non-thermal electrons responsible for the synchrotron radio 
emission. As already noted, an improved theoretical basis 
for the modelling of the radio emission from such systems 
was d eveloped bv lDoughertv et alj (120031) and lPittard et all 
teOOrfl . and we use this as the foundation for our calculations 
of the non-thermal X-ray and 7-ray emission. Key benefits 
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over previous work include a more realistic description of the 
geometry of the WCR, the evolution of the energy distribu- 
tion of the non-thermal electrons as they advect downstream 
from their acceleration site, and absorption of high energy 
photons by pair production in the stellar radiation fields. 
Our predictions for the expected 7-ray emission are based 
on model fits to the available radio and X-ray data, which 
allow us to constrain the range of free parameters and build 
in a degree of consistency for our estimates. 

3.1 IC emission 

IC cooling is a major energy loss mechanism for non-thermal 
electrons in CWB systems, and the non-thermal X-ray and 
7-ray luminosity should exceed the synchrotron radio lu- 
minosity by several orders of magnitude. In the Thomson 
regime, the average energy of a photon after isotropic IC 
scattering is felumenthal fc Gouldlll97(i 

4 2 

hvic ~ -7 hv*. (4) 

Since the average photon energy of an early-type star, 
hv* ~ 10 eV, Lorentz factors of order 10 — 10 4 are suffi- 
cient to produce IC X-ray and 7-ray radiation up to GeV 
energies. The resulting emission has a spectral shape which 
is identical to the intrinsic synchrotron emission at radio 
frequencies. For our calculations, we assume that the scat- 
tering is isotropic and takes place in the Thomson limit. 
While a more thorough treatment would take account of 
Klein-Nishina effects, which are important at high energies, 
recent work has shown that the resulting change to the non- 
thermal electron energy distributio n is rather small when the 
stella r separation is > 2 x 10 14 cm iReimer. Pohl fc Reimerl 
l2006f) . This condition is satisfied in WR140 between orbital 
phases <f> = 0.1 — 0.9, and also in wider binaries like WR 146 
and WR 147. 

Consideration should also be given to the anisotropic 
nature of the IC process, where the emitted power is depen- 
dent on the scattering angle. However, if the O star is in front 
of the WR star (which is the case for WR 140 at cj> = 0.837 - 
see Sec.|lJ, the variation in the resulting emiss ion is less than 
a factor of 2.5 llReimer. Pohl fc Reimer]l2006l) . In Sec. l4~Hl 
we show that the other uncertainties affecting the IC emis- 
sion, such as the value of f e determined from fitting the radio 
data, are at least comparable in magnitude, so we leave this 
issue to future work. 



3.2 Relativistic Bremsstrahlung 

Bremsstrahlung radiation is produced when a charged par- 
ticle is accelerated in the Coulomb field of another charged 
particle. It is possible to obtain emission at 7-ray energies if 
one of the particles is relativistic, since photons of compara- 
ble energy to that of the emitting particle can be produced. 
Inverse bremsstrahlung, produced by collisions between high 
energy ions and thermal electrons, is an insignificant contri- 
bution to the total emission unless the ratio of the energy 
density of non-thermal ions to electrons is very high. There- 
fore, we concentrate on emission produced by collisions be- 
tween non-thermal electrons and thermal elect rons and ions. 
A mo re detailed discussion can be found in iBaring et alJ 
(1999). 



The rate of photon production in the energy interval 
between e 7 and e 7 + de 7 by an electron of kinetic energy E c 



dn~,(E c , e 7 ) 
It 



n e a ee (E c , e 7 ) + > mZfa ep (E e 



(5) 



where e 7 is the photon energy in units of m e c 2 , v e is the 
relative velocity of the colliding particles, and n c are 
the number density of ions and electrons respectively, and 
Zi is the charge of the ion. The electron-proton cross sec- 
tion, (Jep, is the Bethe-Heitler cross section evaluated in the 
Born approximation. We determine a cp for any electron en- 
ergy, relativistic or n on-relativistic, using equation 3BN in 
iKoch fc Motzl Jl959T> . and omit the Coulomb corrections at 
low energies. The cross section for electron collisions with 
fully ionized ions has a Z 2 charge dependence, accounting 
for the factor in Eq.0 For t he electron-electron cross section, 
cr cc , the approach noted in lBaring et alJ lll999F) is followed. 

For 7c 3> 1 and e 7 3> 1, the cross sections depend only 
on the particle's charge and not its mass, so a cp = o ce . The 
ratio of the electron-electron to electron-ion contributions is 
then simply 

(6) 



For material with solar abundances this ratio is « 0.86, while 
in the shocked WR wind of WR 140 it is « 0.47 (see Sec.0. 
At high energies, the bremsstrahlung spectrum (in units of 
photons s _1 cm -2 ) has a spectral index which reflects that 
of the non-thermal particles i.e. p. 



3.3 7T° decay 

High energy hadronic collisions can create neutral pions 
which subsequently decay into two 7-ray photons, e.g., 
p + p — > 7T° + X, 7T — > 7 + 7. The final stage has a branching 
ratio of ~ 98.8% (the decay ir — > 7 + e + + e~ accounts for 
the other 1.2%). The pion decay process yields information 
on the population of non-thermal nucleons, in contrast to 
the IC and bremsstrahlung processes where the emission is 
dependent on the population of non-thermal electrons. 

The cross-section for hadronic collisions depends on the 
ionic species involved, so we first consider the tv° emissiv- 
ity from proton-proton collisions. The production rate for 
neutral pions of energy E^ per thermal proton can be easily 
evaluated using a ^- functional approximation fo r the differ- 
ential cross-section iAharonian fc AtovanllT996Tl : 



8(E^ - f„E k )a pp2 (E)J(E)dE, 



(7) 



where E k is the kinetic energy of the non-thermal protons, 
f„ is the fraction of the kinetic energy of the non-thermal 
proton which is transferre d to the 7r° meson, a PP 2 i s the 
total inelastic cross section jAharonian fc Atovanll996Tl , and 
J{E) is the flux of non-thermal protons of total energy E 
(protons s _1 cm -2 GeV -1 ). Here, only the contribution to 
the 7-ray emission from the neutral pion is considered. After 
the interaction the non-thermal proton keeps approximately 
half of its energy, while the other half goes to pions. On 
average the neutral pion gets approximately one third of 
this energy, or approximately one sixth of the energy of the 
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non-thermal proton, so f n — 0.17. Thus, Eq. |7| becomes 
(pionss^GeV^H" 1 ): 



qn(E^) = —a PP 2(E)J(E). 
J™ 

The 7-ray emissivity (photons s _1 GeV" 1 !!" 1 ) is then 



q 1 (E 1 



VeJ- 



- r dE„, 



(8) 



(9) 



where i5 m i n = E^ + m^ c /4£ 7 and m 



is the mass of the 

neutal pion (134.9766 MeV/c 2 ). 

In CWB systems, we must account for the composition 
of non-thermal and thermal ions. For collisions involving two 
nuclei with atomic weigh ts Ai and A2, the cross section is 
JOrth fe Buffingtonlll97g) 



3 / 8 + ^3/8 



1) 



CTpp2. 



(10) 



This prescription is appropriate for proton-helium interac- 
tions, but its validity for heavy nuclei, such as Fe, is un- 
clear. It is adopted here since its potential inaccuracy is 
small compared to the other uncertainties and assumptions 
in this work. In the case of WR stars, such as WR 140, the 
wind contains little or no hydrogen, and the cross-section 
for pion decay is best calculated relative to the He-He cross 
section a aa = 5.6<t PP 2. 

Considering only protons and Helium nuclei, an en- 
hancement to the pion production rate in Eq. |H] of 1.61 is 
obtained for material of solar composition. Greater enhance- 
ments are obtained for the nuclear processed material in WR 
winds. For example, if collisions involving He, C, O and Ne 
are considered for the WC7 wind abundances in WR 140 
(see Sec. l4.2.Tl . then a ~ lAa aa ~ 7.8cr pp 2. 



3.4 7-ray absorption 

7-rays above a suitable threshold energy may create 
electron-positron pairs via interaction with a lower energy 
photon or with charged nuclei. While the latter process has 
an insignificant optical depth in CWB systems, in contrast 
the high stellar radiation energy density may be a significant 
source of opacity to high energy 7-rays. The probability of 
absorption depends on the cosine of the angle between the 
directions of the two photons, fi. For a high energy photon 
of energy E-y interacting with a stellar photon of energy e, 
the optical depth is given by 



dr 
dr 



°"(x)(l - ^)n(e, H, r)ded(J,, 



(11) 



where x = \J -E 7 e(l — n)/2 is the centre-of-momentum 
frame energy 2 scaled by m e c 2 , 



1 i nr <= 
X 6 



(2 X 4 + 2 X 2 - 1) ln( X + Vx 2 - 1) 



-xa+x^vx 2 



(12) 



where r c = e 2 /r 



and n(e, [1, r) is the number density of 



stellar photons. The threshold for pair production, and thus 
for absorption, is \ ^ 1. 

The equation for \ contains a typo in iBaring fe Hardind tl99^ . 




Figure 2. The optical depth from two-photon pair production 
as a function of the non-thermal photon energy for a single line- 
of-sight from the stagnation point of the WCR in WR 140. The 
shape of the curves reflects the Planck distribution for the stellar 
photons (T = 45,000 K). We show calculations at 4> = 0.837 for 
models A and B (see Sec.|3J. In each case, only the photons from 
the O star are considered. The energy of peak absorption is a 
function of the viewing angle, 8, and increases as 8 decreases. 



In wide CWBs it can be assumed that the stars are 
point-like, and that the stellar photons stream radially away 
from each star. This allows the (1 — /i) term to be taken out 
of the integral in Eq. 1111 The optical depth as a function 
of non-thermal photon energy from the stagnation point to 
an observer for two models of WR 140 is shown in Fig. 
Table summarizes some key parameters of these models. 

Of the two models shown in Fig. [5] both the shock apex 
and the line of sight from the apex are closer to the O star in 
model B. The greatest opacity occurs at ~ 100 GeV, where 
optical depths as high as ~ 25 (model B) may be reached. 
At lower and higher energies the optical depth decreases, 
and optical depth unity occurs at approximately 50 GeV 
and 5 TeV for model A, and at approximately 20 GeV and 
40 TeV for model B. In practice, the effective optical depth 
is likely to be somewhat smaller due to the spatial extension 
of the emission region. More accurate calculations using ra- 
diative transfer along multiple sight lines are presented in 
Sec. gl] The CWB system WR 146, which has a wider stel- 
lar separation, is expected to have a maximum optical depth 
to pair-production of ~ 1. In WR147, which is wider still, 
the attenuatio n is negligible. Similar conclus ions have also 
been noted bv lReimer. Pohl fc Reimerl (1200611 . 



4 NON-THERMAL EMISSION MODELS OF 
WR140 

In this section we apply our model to WR 140 (HD 193793), 
the archetype of long-period CWB systems. WR 140 con- 
sists of a WC7 star and an 04-5 star in a highly el- 
liptical orbit (e ~ 0.88), and is well known for the 
dramatic variation s in its emission fro m near-IR to ra- 
dio wavelengths iWilliams et alJ Il990t IWhite 



Beckerl 

11995 ), and also at X-ray energies during its 7.9-year or- 
bit Zhekov fc Skinned l2000l: iPollock. Corcoran fc Stevens! 
l2002t IPollock et al.ll2005h . The variability appears to be 
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linked to the WCR, which experiences significant changes as 
the stellar separation varies between ~ 2 AU at periastron 
and ~ 30 AU at apastron. The observed radio emission in- 
creases by up to two orders of magnitude between periastron 
and a frequency-dependent peak between orbital phases 0.65 
to 0.85, followed by a steep decline. The decline of the radio 
emission prior to periastron is consistent with the extinc- 
tion suffered by the X-ray emission at these orbital phases. 
The X-ray lightcurve does not follow the expected 1/D vari- 
ation for an adiabatic WCR dPollock. Corcoran fe Stevens! 
120021) . which may i ndicate that the WCR becomes ra dia- 
tive near periastron dVarricatt. Williams fc Ashokll2004T) . It 
is also possible that this discrepancy is caused by a sub- 
stantial amount of the available energy in the wind-wind 
collision being placed into non-thermal particles, with the 
efficiency varying as a function of phase, though we believe 
this scenario should be considered as less likely. 

WR 140 is also noteable for its possible association 
with an unidentified EGRET source, lying on the out- 
skirts of the positional error bo x of 3EG J2022+4317 
iRomero. Benaglia fc Torres] ll99ST) . The error box has a 
radius of w 0.7°, but WR140 is the only known high- 
energy source in the vicinity. 3EG J2022+4317 has a flux of 
(24.7±5.2) x 10 -8 photons s -1 cm -2 , and a photon spectral 
index V = 2.31±0.19 (where N(E) oc £T r ) Martman et~ail 
Il999h . The va riability of 3EG J2022+43 17 is discussed in 
some detail bv lBenaglia fc Romero! a2003l) . Since it was only 
detected in 3 out of 12 observing periods (at approximate 
orbital phases 0.97, 1.00, and 1.16), the upper limits of the 
remaining observations may hide real variability. 

It has been suggested previously that better agree- 
ment may be obtained between observed and modelled ra- 
dio lightcurves if the wind of the WR star is disc-like 
JWhite fc Beckerlll995h . but a specific (and unusual) orien- 
tation is re quired for consistency w ith polarization measure- 
ments (see iMarchenko et aljl2003l and discussion therein), 
and t here is no evidence for such a disk in X-ray observa- 
tions JPollock. Corcoran fc StevenslEio^ . In the following, 
it is assumed that the undisturbed wind is isotropic. 

4.1 Details of the model 

Our calculations of the observed thermal and non-thermal 
emission from CWB systems are based on the density and 
temperature distribution obtained from 2D hydrodynami- 
cal simulations of the s tellar winds and th eir collision (see 
iDouehertv et alJ l2003t IPittard et afll2006ft . At radio fre- 
quencies, both the stellar winds and the WCR contribute 
thermal emission and absorption. The free-free emission and 
absorption coefficients are determined from the local tem- 
perature and density values on the hydrodynamic grid, and 
clumping is accounted for (see Sec. 14.2.41 . The non-thermal 
emission is modelled after specifying the population and 
spatial distribution of non-thermal particles, as noted in 
Sec. [5] Additional details relating to the hydrodynamical 
code, and the method and assumptions used to obtain the 
various emi ssion and abso r ption coefficients in each cell, can 
be found in lPittard et alJ i2006T) and references therein. 

The assumption of axisymmetry in our 2D mod- 
els is good until orbital phases of ~ 0.97, as 
the orbi tal eccentricity is so high. S ince radiative 
braking dCavlev. Owocki fc Cranmerl 11997) or inhibition 



JStevens fc Pollocklll994T) is not expected to be important, 
the wind speeds are assumed to be constant. The calcula- 
tions are performed on a grid with 1600 x 800 cells, which 
spans the range < r < 4 x 10 15 cm, —4 x 10 15 < 2 < 
4 x 10 15 cm, and is large enough to capture the vast major- 
ity of the radio and X-ray emission from the WCR, though 
inevitably the models underestimate the radio flux at the 
lowest frequencies considered. The line of sight into the sys- 
tem is specified by the angle 9, which is measured from the 
line perpendicular to the line of centers of the two stars: 
9 = +90° corresponds to the O star in front of the WR 
star, 9 = 0° corresponds to quadrature, and 9 = —90° cor- 
responds to the WR star in front of the O star. 



4.2 Parameters of WR 140 

Recent high-resolution VLBA observations at 8.4 GHz be- 
tween orbital phase 0.7 and 0.9 have provided important new 
constraints to models of WR140 iDoughertv et al.ll2005T) . 
An arc of emission is observed, resembling the bow-shaped 
morphology expected for the WCR, and rotates as the or- 
bit progresses. This rotation allows derivation of the orbit 
inclination and semi-major axis, leading to a geometric dis- 
tance determination (1.85T0.16 kpc) that is independent of 
stellar parameters. This distance is greater than previously 
thought, and implies the O star is a supergiant. Therefore, 
we adopt log (Lboi/Lo) = 6.18 and Woo = 3100 km s -1 
for the 04-51 primary, and log (Lboi/ Lp) = 5.5 and 
«oo = 2860 kms -1 for the WC7 secondary dDoughertv et afl 
I2005T) . These lumin osities are higher than averag e for 
the spectral types. liTerrero. Puis fc Naiarrol i2002T) and 
iRepolust et alJ <|2004F) suggest a typical value of 5.7-5.9 for 
an 04-5 supergiant, while the mean luminosity for WC7 
stars is about 0.4 dex lower (although the scatter in the ab- 
solute magnitude of WC7's is ~ 1 magnitude). Therefore, 
there may be some scope for lowering the adopted luminosi- 
ties. 

[Dougherty et all 1)20051) used the models of 
IZhekov &H5khmBrl 11200(1) of the observed thermal X- 
ray luminosity to derive a mass-loss rate for the WR 
star of 4.3 x 10 -5 M©yr -1 . However, there are several 
complicating factors, such as the inadequate constraints 
on the composition of the WC7 wind and on the wind 
momentum ratio, rj. As n increases, the opening angle of 
the WCR increases, and a greater fraction of the WC7 wind 
is shocked, resulting in an increase in the X-ray luminosity 
if the mass-loss rate of the WR star is unchanged. There 
is also the possibility of some degree of shock modification. 
These issues are explored in the following sub-sections, 
leading to improved estimates of the mass-loss rates and r\ 
in WR 140. 



4-2.1 WR wind composition 

The ther mal X-ray emission is larg ely from the shocked WR 
star wind jPittard fc Stevensl2002r) . Hence, it is particularly 
sensitive to the abundance of He and C since the free elec- 
trons are largely stripped from these elements. For a pure He 
and C wind, the continuum flux scales as (2 2 + (C/He)6 2 ), 
whe re C/He is the rati o of abundance by number of C and 
He llPollock et al .120051) . Unfortunately, the mass fraction of 
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Table 1. Some key parameters for our models of WR140 at 
<j> = 0.837 (Sec. |4j. D is the binary separation, r\ is the wind 
momentum ratio, tq is the distance of the stagnation point of 
the WCR from the O star in units of the stellar separation D, 
and d is the distance of the system. [/ t hi Ub, U p h, n e , and B are 
all evaluated at the stagnation point, while the pre-shock Alfvcnic 
Mach number, Ma, is evaluated immediately upstream of the WR 
shock on the line of centers, assuming that at this point the shock 
is highly perpendicular and that the B-field is reduced by a factor 
of 4 from its post-shock value. 8 is the angle of the line of sight 
into the system. Two values are given for n e - these are for the 
WR and O star side of the contact discontinuity respectively. 



Parameter 


Model A 


Model B 


D (10 14 cm) 




3.13 


lo S ( L o,boi/L Q ) 




6.18 


log (£ W R,bol/L ) 




5.5 


d (kpc) 




1.85 


n 


0.22 


0.02 


ro 


0.32 


0.124 


Crel,e 


1.38 x 10" 


' 3 5.36 x 10~ 3 


Cb ' 


0.05 


0.05 


U t h (erg cm -3 ) 


0.78 


1.2 


Ub (ergcm~ 3 ) 


0.04 


0.06 


f^ph (erg cm" 3 ) 


1.6 


10.3 


n e (WR) (cm" 3 ) 


3.7 x 10 6 


6.0 x 10 6 


n c (OB) (cm" 3 ) 


1.3 x 10 7 


2.0 x 10 7 


B (mG) 


990 


1200 


M A 


11.4 


11.9 


e 


5° 


35° 



C in WC stars spans the range 0.2 — 0.6, with no strong 
sensitivity with spectral type. In this work, we assume that 
C/He = 0.25 by number (C= 0.4 by mass), as the ratio 
of Civ 5471 /He ii 5411 is comparable to that in WR90 
JPessart et alj2000l. P. Crowth e r, priv . communication). We 
note that lEenens fc Williams! <ll992l) had previously mea- 
sured C/He = 0.15 for WR 140. The O abundance in WC 
stars is also poorly constrained, though the number ratio 
of O/He is esti mated as 0.03 ± 0.01 for W R90 (WC7) and 
WR 1 35 (WC8) l|Dessart et al.ll200ftl . Since IffiUier fc Millerl 
J1999T had earlier obtained O/He=0.1 for WR111 (WC5), 
we adopt O/He = 0.05 (O = 0.1 by mass). Our val- 
ues for_C/He_andO/He lie within the range considered 
by IPollock et alJ (2005), while we note that in comparison 
l^reitov^g^kmnBn(l2OToi) adopt C and O abundances which 
are below the lower end of the ranges noted above. With 
these assumptions, the mass fractions of hydrogen, helium, 
and "metals" (X, Y, and Z respectively) for the WC7 wind 
are X = 0.0, Y — 0.5, Z = 0.5. Solar abundances are as- 
sumed for the O star (mass fractions X = 0.705, Y — 0.275, 
Z = 0.020). Both winds are assumed to have a temperature 
of 10,000 K, with an ionization structure of H + , He + , and 
CN0 2+ . 



4-2.2 Collisionless shocks - electron thermalization, 

non- equilibrium ionization, and shock modification 

Previous models of the X-ray emission from the WCR 
in wide CWB systems have usually assumed that the 
shocks are collisional, with immediate equilibration of the 



post-shock electron and ion temp eratures. However , in re - 
ality they are collisionless (e.g., iDraine fc McKeel Il993l) . 
and there is growing observational evidence that in such 
shocks the ratio of post-shock electron to ion te mperature, 
T e /T j, is a function of the shock speed (e.g., iRakowskil 
120051) . Temperature equilibration through Coulomb colli- 
sions and plasma instabilities then occurs some distance 
do wnstream. This process h as been examined in WR 14 
bv IZhekov fc Skinnerl l)2000h and bv IPollock et alJ l)2005lh 
and can naturally account for the observed soft X-ray con- 
tinuum emission. Since the ionization timescale is of the 
same order as the timescale for energy transfer from ions to 
electrons, it is not surprising that observations of WR 140 
also indicate that the post-shock plasma h as yet to reach 
ionization equilibrium ijPollock et alJl2005l ). and that line- 
profile models which assume rap id ionization equilibrium 
jHenlev. Stevens fc Pittardl l2003f) are unable to reproduce 
the observed c orrelation of X-ray line widths with ioniza- 
tion potential jHenlevll2005l) . 

Despite this progress, the physics of the heating process 
in the shock layer remains poorly understood, and shocks 
which are efficient at accelerating particles may equilibrate 
their electron and ion thermal populations more rapidly, 
in addition to reducing the equilibri um value of T c and T; 
fellison. Decourchelle fc Balled liooi) . In this case, the X- 
ray spectrum is softened due to the reduction in the post- 
shock electron temperature behind a weaker subshock (see 
Fig. 0. It is unclear whether, and to what extent, shock 
modification may occur in WR 140, but this mechanism 
could potentially play some role in creating the low value of 
T e observed. A simple calculation reveals that even if parti- 
cle accleration is highly efficient, the non-thermal particles 
in the precursor would have a negligible effect on the ioniza- 
tion of the pre-shock gas, consistent with the observed lack 
of lin e emission near the wind terminal speeds ijPollock et alJ 
120051) . 

A full understanding of WR 140's spectrum will even- 
tually require a detailed hydrodynamical model which ac- 
counts for the large post-shock equilibration timescales of 
the plasma ionization and electron and ion temperatures, 
combined with the acceleration of non-thermal particles and 
their effect on the shocks bounding th e WCR and both th e 
upstream and downstream flow (e.g., IPollock et alJlioOSl ). 
As the development of such models is not trivial, here the 
synthetic X-ray spectra are calculated from models of un- 
modified collisional shocks. Given this simplification and the 
uncertainties described above, we adopt a rather pragmatic 
approach, and concentrate on obtaining a comparable flux at 
~ 3 keV. Absorption is negligible at this energy at the chosen 
phase (see Fig. I18H for reasonable viewing angles, while at 
higher energies the relative softness of the observed emission 
is clear evidence for T c < 2], with the (less likely) possibility 
that shock modification also plays a role (see Fig.^J. 

The fact that the thermal X-ray luminosity from the 

• 2 

WCR, L x , scales as M when the post-shock gas is largely 
adiabatic means that even if shock modification occurs, and 
were to lead to fairly substantial (e.g., factor of two) changes 
in L x , mass-loss rates derived from models with unmodified 
shocks will be much less affected. Moreover, the inferred 
mass-loss rates are largely insensitive to the assumption of 
T c = Ti. For instance, setting T c = 0.571 throughout the 
entire WCR reduces the emission at 3 keV by 23% com- 
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pared to a model with T c — 71. An increase of only 14% 
in the mass-loss rates is required for the original flux to be 
regained. In this way we can still obtain useful estimates of 
the stellar mass-loss rates from the modelling work which 
follows. 



4-2.3 WR mass-loss rate from thermal X-ray models 



Having specified the composition of the winds (Sec. 14. 2. II 
and unmodified, collisional shocks (Sec. 14. 2. 2^ . theoreti- 
cal calculations of the X-ray emission, based on axisym- 
metric 2D simulations of the stellar winds and the WCR, 
can now be made. The X-r ay emissivity is calculated using 
the MEKAL emission code (iMewe. Kaastra fc Liedahll995l 
and references therein) for optically thin thermal plasma, 
and ionization equilibrium and rapid thermalization of the 
electrons is assumed. 

The intrinsic (i.e. prior to any absorption) X-ray emis- 
sion calculated from a model of the WCR with a fixed WC7 
mass-loss rate, but for two different values of rj, is shown in 
Fig- El As r\ increases, the X-ray flux increases, as explained 
earlier. Clearly, in order that the X-ray flux calculated from 
hydrodynamical models of the WCR is comparable to obser- 
vations, the mass-loss rate of the WC7 star must be adjusted 
as r\ is varied. 

There have been a number of X-ray observations 
of WR 140 in recent years, and we choose to use an 
ASCA observation taken on Dec. 20 1999 (<j> = 0.837). 
At this phase the X-ray luminos ity is climbing rapidly 
JPollock. Corcoran fc Stevens! I2002T) . while the radio emis- 
sion is around its maximum intens ity, and the stellar sep - 
aration, (D = 3.13 x 10 14 cm, see IDoughertv et al.l 120051) . 
is sufficiently large that the shocked gas in the WCR is 
approximately adiabatic, allowing easy adjustment of the 
model mass-loss rates in order to obtain a specific X-ray 
flux. The mass-loss rates required to match the observed X- 
ray emission at <f> = 0.837, as a function of rj, are noted in 
Table [2] Fig. [|] demonstrates that the X-ray emission from 
two of these models is comparable to the observed emis- 
sion at ~ 3 keV, as required. Models of the radio emission 
(Sec. l4.3"^t indicate r\ is likely to have a value near the lower 
end of the range covered by Table [~] This implies the mass- 
loss rate of the 04-5 supergiant in WR 140 is lower (~ 8 x 
10~ 7 Mp) yr _1 ) than the suggested values in lRepolust et all 
J2004h (8.6 - 8.8 x 10~ 6 Meyr -1 ). However, an order of 
magnitude reduction i n the mass-loss rate is consistent w ith 
the recent estimat es oflFull erton. Massa fc Prinial i2006f) . 

We note that iPollock et alj (120051) determine M W r ~ 
1.5 x 10 -5 M© yr _1 from the degree of circumstellar absorp- 
tion observed by Chandra near periastron. While this value 
is independent of the clumping in the winds, it is likely to be 
a lower limit since the overall absorption of emission from 
the WCR is probably somewhat less than occurs for emission 
near the apex of the WCR. This method will also underesti- 
mate Mwr if there are regions of O star wind between the 
observer and the WR star at this phase due to the curva- 
ture of the WCR induced by the orbital motion of the stars. 
Therefore, it is argued that the range of values of Mwr con- 
sidered in our models is consistent with this constraint (see 
Table ["}, and does not preclude values of rj as low as 0.02. 

Due to the relatively low abu ndance of C and O as- 
sumed bv lZhekov fc Skinned ll200(J) . their derived mass-loss 




Figure 3. Comparison of the intrinsic thermal X-ray emission 
calculated from hydrodynamical models of the WCR in WR 140 
at <j> = 0.837 for two different values of rj, assuming unmodified 
collisional shocks (Sec. 14.2.21 . r? = 0.22 (top) and rj = 0.0353 
(bottom). In these calculations the mass-loss rate of the WC7 
star and the terminal wind speeds were fixed (Mwr = 4.3 x 
10- 5 M Q yr- 1 , Docwr = 2860 kms" 1 , i^q = 3100 kms" 1 ), 
while the mass-loss rate of the O star was varied according to 
the desired value of r\. The continuum flux is also sensitive to the 
adopted abundance of the WR wind (see Sec. I4.2.T1 . 



Table 2. The mass-loss rates (in units of 10 — 6 Mq yr -1 ) required 
to match the observed X-ray flux at 3 keV and <j> = 0.837 as a 
function of r). Column 4 indicates lower limits for the volume 
filling factor, /, of clumps in the winds, as determined from a 
comparison between the X-ray and radio derived mass-loss rates 
fSec. l4.2~4t . For simplicity, we assume that the two winds have the 
same values of /. Recent analyses of O star wind s indicate that 
/ ~ 0.02-0.1 fsee lFullerton. Massa fc Primateo"""". and references 
therein). Columns 5-7 note approximate values of 9 for which the 
sight lines are parallel to the asymptotic opening angles of the 
WR shock, contact discontinuity, and O shock. 
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rates for rj = 0.0353 are higher than ours (cf. their model A 
and AA), even though the system distance has since been 
revised upwards. We also note that the wind parameters 
noted in IDoughertv et all l)2005ll overpredict the observed 
X-ray emission by a factor of about 8. 

Given the existing uncertainties with respect to the de- 
gree of shock modification, and the amount of post-shock 
electron heating and the timescale for electron and ion tem- 
peratures to reach equilibrium, Fig. [|] demonstrates that it 
is difficult to uniquely determine the WR star mass-loss rate 
and r\ from the X-ray spectrum alone. While a simpler, and 
perhaps a more reliable, method of determining r\ in CWB 
systems with highly eccentric orbits involves utilizing the du- 



12 J. M. Pittard and S. M. Dougherty 



1e-08 



1e-09 



1e-12 





ASCAcHO.837 




Intrinsic o=0.22 


j 

i 

j 

.1 

ill. * 

■:■ Jj ... 


Intrinsic n=0.02 

; i 

!!■.:,] 

!•■■!• a. . 









i I 1 1 1 — I — I — I — ^-1 1 . 1 — I — I — I — 

0.1 1 10 

Energy (keV) 



Figure 4. The intrinsic thermal X-ray emission from hydrody- 
namical models of the WCR at orbital phase <j> = 0.837 with 
n = 0.22 (dashed) and n = 0.02 (dotted). There is a slight de- 
crease in the hardness of the model spectra with n, as expected. 
Also shown is the observed X-ray emission (solid). The X-ray 
emission from the models (at 3 keV) is comparable to the observed 
emission, though the latter is softer. As discussed in Sec. 14.2.21 
this is primarily because T e < 7]. 



ratio n of any X-ray "eclipse" (see, e.g.. IPittard fc Corcoran! 
| 2002|) , the very shor t eclipse seen in the lightcurve of WR 140 
JPollock et alj2005Tl means that a detailed analysis using 3D 
hydrodynamical simulations is probably needed. Therefore, 
we turn to the high quality radio spectra now available to 
further constrain 77 (Sec. l4.3"!?l . 

4-2.4 Wind clumping 

There is now considerable observational e vidence that the 
wind s of massive stars are clumpy (e.g., |Crowther_et_^ 
20021; iHillier et al.ll2003l: iMassa et al.1 120031: [Repolust et al, 
2004; iFullerton. Massa fc Prinial 120061 : lOwocki fc Cohen 
2006). Additional support is provided by models of 
radiatively-driven winds. Clumping affects the thermal ra- 
dio emission from the unshocked winds, though simple ar- 
guments suggest that in wide binaries, where the WCR is 
largely adiabatic, clu mps are rapidly d estroyed upon en- 
countering the WCR JPittard et alJl200d) . Therefore, in the 
following models we assume that the unshocked winds are 
clumpy (with the same degree of clumping for both winds), 
and the WCR smooth. 

The degree of clumping may be conveniently parame- 
terized by a volume filling factor, /, where it is assumed 
that the clumps occupy a fraction / of the wind volume, 
with the interclump medium being of negligible density. 
Radio-derived mass-loss rates of clumpy winds are related 
to those of smooth winds by M = M smooth// 1 ' 2 , since the 
thermal flux Ss in creases accordingly (see, e.g., Eq. 16 in 
IPittard et al .120061) . This relationship is more complicated in 
binaries, as the degree of clumping may be different for each 
wind, and the thermal emission from the WCR may make 
a non-negligibl e contribution to the total thermal emission 
(see Fig. 11 in IPittard et al.ll2006ft . Nevertheless, we pro- 
ceed by assuming the WR wind dominates the thermal radio 
emission in WR 140, so that a value for / can be deduced 
from a comparison of the X-ray and radio derived mass-loss 



rates. The unclumped estimate of 5.3 x 10 5 Mq yr 1 from 
IWilliams et al] Jl99Cj) translates into 9.7 x 10" 5 M yr -1 for 
a system distance of 1.85 kpc, and leads to the lower limits 
for / noted in Ta bled An additional un certainty is that the 
mass- loss rate of IWilliams et al] (Il990l) assumes the radio 
flux at <f> = 0.0 is entirely thermal, which is far from clear. 
While such assumptions should be avoided if at all possible, 
our approach is conservative in the sense that the theoretical 
thermal flux from the unshocked winds cannot exceed the 
measured flux at periastron. Since the total observed flux at 
4> = 0.837 is ~ 20x greater at 5 GHz and ~ 5x greater at 
22 GHz than our estimate of the thermal flux, the models in 
Sec. 14.3.21 should not be overly sensitive to our choice of /, 
although p will be overestimated slightly if the thermal flux 
is in fact lower than assumed. Wind clumpi ng also affects 
the r adius of the r = 1 surface (see Eq. 17 in IPittard et al] 
120061) . In the clumpy WR wind in our models, J? T =i has 
values of 8.8, 4.0, and 1.8 x 10 14 cm at 1.6, 5, and 15 GHz. 
This is independent of the variation in Mwr with r\ due to 
the corresponding variation in /. 

4.3 Models of the radio data 

We concentrate on obtaining reasonable spectral fits to the 
radio data at the same orbital phase (<f) = 0.837) we deter- 
mined potential values for stellar mass-loss rates and corre- 
sponding values of 77 from the observed X-ray emission (see 
Sec. 14. 2. 3"1 . The brightness temperature of the radio emission 
is ~ 3 x 10 7 K at this phase, so the o bserved emission i s dom - 
inated by non-thermal processes tooughertv et al.l 12005^ 1. 
Therefore, a good estimate of the population of relativis- 
tic electrons can be obtained, allowing predictions of the 
non-thermal X-ray and 7-ray emission. We will apply our 
modelling to a variety of orbital phases in a forthcoming 
paper. 

There are several free parameters used to fit the ob- 
served radio emission, including 77, the index of the non- 
thermal electron energy distribution, p, and £ re i, e and £b, 
which normalize the model spectra and the magnetic field 
strength (and hence the impact of the Razin effect). The 
effect of reducing the luminosity of the stars is also investi- 
gated. The line of sight angle into t he system at (f> = 0.837 is 
specified by the orbit derivation of iDouehertv et al] ]2005l) 
as 8 — 52°, though the uncertainty in this value is large 
(~ 10°) and we allow 9 to vary substantially (±10 — 20°) in 
our models. 



4.3.1 Parameter study 

Before attempting a fit to the (f> — 0.837 data, it is instruc- 
tive to examine how the model spectra vary with 8, 77, p, £b, 
and Lboi- Fig. [5] shows how the model synchrotron spectra 
change as 9 (and hence the free-free absorption along the 
line of sight) is varied, for models with r\ — 0.11, p — 2.0, 
Crei.e = 0.215, and Cb = 2.6 x 10~ 4 . The Razin effect pro- 
duces a turnover at about 2.3 GHz, while the high fre- 
quency emission from each model tends to the same slope 
as the free-free attenuation becomes negligible. The sensi- 
tivity of the spectrum at low frequencies to changes in 8 
is expected since the size of the radio photosphere is larger 
than the stellar separation. At high frequencies the radio 
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photosphere is smaller and has a commensurately weaker 
influence on the absorption of emission from the WCR. The 
synchrotron emission is least affected by free-free absorption 
when 9 m 20° (i.e. when the lines of sight are roughly parallel 
to the asymptotic angle of the WR shock) , and the turnover 
frequency is the same as that of the intrinsic synchrotron 
emission. A dramatic increase in the level of free-free ab- 
sorption occurs as 9 decreases and the WR wind moves in 
front of the WCR, with maximum absorption at 9 = —90° 
(i.e. WR star in front). For 9 < 20°, free-free absorption 
causes the low frequency turnover. As the line of sight swings 
through the WCR from the WR shock to the O shock i.e. 
9 — 20 — > 56°, the level of absorption again increases, al- 
though in a somewhat complicated fashion as the source is 
extended and there are multiple lines of sight. While free-free 
absorption determines the slope of the spectrum between 3 
and 20 — 30 GHz, it is only between 50° and 60°, when 
the sight lines move out of the WCR and into the O wind, 
that the turnover due to free- free absorption dominates over 
the Razin effect, and the turnover becomes very broad. The 
turnover frequency increases as the O wind moves in front 
of the WCR i.e. as 9 — > 90°. The smaller stellar separation 
in the models presented here results in free-free absorption 
having an influence to muc h higher frequenc i es tha n in the 
models shown in Fig. 10 of iDoughertv et al.l (120031) . 

It is interesting to note that when the sight lines are 
through the WCR, a turndown or "kink" is visible in the 
spectrum at a few 10's of GHz (this is most clearly seen 
when 9 — 50°). Typically, such turndowns are attributed to 
a change in the non-thermal energy spectrum. However, it 
is certainly not the cause here, because the intrinsic syn- 
chrotron spectrum does not display such a kink (although 
spectral breaks from IC cooling occur in our models, the 
break frequency is spatially dependent, and is smoothed out 
o nce the flux is integr ated over the entire WCR - see Fig. 3 
in iPittard et aTll200rJ) . It is clear that lines of sight through 
the WCR are optically thin at frequencies above 20-30 GHz, 
and almost all, if not all, the synchrotron emission from the 
WCR is seen. 

The variation of the spectra as a function of r\ is shown 
in Fig. HJ where the total emission (synchrotron plus free- 
free) is displayed. In this and the following figures, emission 
between 0.5 and 50 GHz is examined, since this is approxi- 
mately the frequency range of the observations. The spectral 
variations in Fig. |5] bear some resemblence to the changes 
with 9, since as r\ varies, the relative angle of the sight lines to 
the WCR changes if 9 is fixed. For ? 7 = 0.055 - 0.16, the line 
of sight is through the lower opacity WCR. When r\ = 0.22, 
it is through the unshocked O wind, and the synchrotron 
emission from the WCR suffers more free-free absorption 
than seen in the other models, which results in the low fre- 
quency turnover occuring at a higher frequency (c.f. Fig. |5J. 
For the remaining models, the reduced level of free-free ab- 
sorption means that the turnover is the result of the Razin 
effect, and therefore occurs at approximately the same fre- 
quency in each case (since £b is fixed). It is clear from Fig. El 
that the spectral slope between 10-40 GHz steepens as r\ de- 
creases (the slight upward curvature in the rj — 0.055 model 
near 40 GHz is caused by rising free-free emission), which 
is due to the change in the relative angle of the sight lines 
with respect to the WCR, as seen in Fig.0 and the resulting 
change in the free- free opacity. While Fig. |S| does not show 



frequencies as high as Fig. [S] the synchrotron emission tends 
to a similar slope in all models (since p is the same in each 
model) with a slight steepening as rj decreases due to the 
WCR moving closer towards the O star, and the amount of 
IC cooling increasing. 

Fig. |7| shows the effect of varying the slope of the non- 
thermal electron energy distribution. As p decreases, the 
slope of the intrinsic synchrotron emission increases, result- 
ing in flatter attenuated spectra. Fig. |5] illustrates the ef- 
fect of changing £b- As £b is reduced, the spectral index 
above 4 GHz increases. By decreasing £b (and thus B), the 
synchrotron emission at a given frequency arises from non- 
thermal electrons with higher Lorentz factors (y oc 7 2 -B), 
and thus sample a steeper slope for the non-thermal electron 
energy distribution due to the increased impa ct of IC cool- 
ing at high 7 (see Fig. 3 in lPittard et al.ll2006t) . At very low 
values of £b, the Razin effect becomes the dominant mech- 
anism responsible for the low frequency turnover in Fig. |H] 

The effect of adjusting the stellar luminosities is shown 
in Fig. |5| For simplicity, the luminosity of both stars is ad- 
justed by the same amount, although the amount of IC cool- 
ing is primarily determined by the luminosity of the O super- 
giant. In our model, a power-law distribution of non-thermal 
electrons is injected along the entire shock front with the 
same initial spectral index and normalization (relative to 
the local post-shock thermal energy density - see Sec. 12. 2, It . 
Near the apex of the WCR, there is strong IC cooling but 
the cooling weakens away from the apex as the stellar radia- 
tion field density diminishes, and the non-thermal particles 
can maintain higher energies. As the Lboi values are reduced, 
the influence of IC cooling decreases, the non-thermal elec- 
trons are able to keep higher energies throughout the WCR, 
and the population of the highest energy electrons can exist 
throughout a larger part of the WCR. The result is that the 
spectrum of synchrotron emission, summed over the entire 
WCR, becomes less steep. Clearly a specific spectral slope 
can be obtained with a larger value of p if the model values of 
Lboi are reduced. For a given energy density of non-thermal 
electrons immediately post-shock, a reduction in IC cooling 
also produces an increase in the overall flux. It is conceiv- 
able that the assumptions in Sec. I4.2l mav overestimate the 
stellar luminosities, but with one exception the models in 
the rest of this work are based on the luminosities noted in 
Table 

We have also investigated what happens to the observed 
synchrotron emission if the emission region is limited to r < 
7rro/2, as suggested bv lEichler fc Usovllll993l) . We find that 
the spectral index is essentially unchanged, and thus our fits 
for p are robust in this context, though the flux is diminished 
to 40 — 75% of its previous value due to the reduced size of 
the emitting region. 

4.3.2 Spectral fits at <j> = 0.837 

Our inital aim was to determine parameter values directly 
from model fits to the radio data. Unfortunately, reasonable 
fits to the data are possible with a range of parameters. For 
example, by varying p and 9 similar fits are attained over a 
wide range of rj. Furthermore, the data can be fit with either 
free-free absorption or the Razin effect responsible for the 
low frequency turnover. However, despite this degeneracy 
of models, it is possible to discriminate between models by 
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Figure 5. The influence of 9 on the observed synchrotron emis- 
sion. All models have r\ = 0.11, p = 2.0, Crel.e = 0.215, and 
Cb = 2.6 X 10~ 4 . The lines of sight are through the WR wind 
when -90° < < 19°, through the WCR when 19° < 6 < 56°, 
and through the O wind when 56° < < 90°, as indicated by 
the long-dashed, short-dashed, and dotted lines respectively (see 
Table |2J. The intrinsic synchrotron emission is indicated by the 
solid line. 



Figure 6. The influence of rj on the model radio spectra. IC 
cooling and the Razin effet are included in these calculations. 
Ionic cooling and synchrotron self-absorption are negligible. All 
models have = 50°, p = 2.0, Crel.e = 0.215, and Cb = 2.6 X 
10 -4 , though the mass-loss rates and / were varied according 
to Table |2] The line of sight is through the WCR in all models 
except for rj = 0.22, when it is through the unshocked O star wind. 
The synchrotron emission from the WCR dominates the free-free 
emission from the winds and WCR for the chosen parameters. 



using other properties, such as the values of 6 and C re i,i re- 
quired by each model, which we discuss later in this section. 
Parameters for a variety of model fits are listed in Table 03 
Two model fits which span the range of r\ considered 
are shown in Fig. 1101 For both of the models in Fig. 1101 
Cb = 0.05 was assumed. Sightlines which have 9 marginally 
less than #wr are necessary for free-free absorption to cause 
the low frequency turnover. With these restrictions, p — 1.4 
is needed to match the slope of the data points between 5 
and 22 GHz. Fits with free-free absorption dominant over 
Razin cannot be obtained when the sightlines are through 
the WCR because the free-free absorption is too low, while 



Figure 7. The influence of p on the model radio spectra. All 
models have r, = 0.11, = 50°, Crel.e = 0.215, and Cb = 2.6 X 
10 — 4 , and the line of sight is through the WCR. 



Figure 8. The influence of Cb on the model radio spectra. All 
models have rj = 0.0353, p = 1.4, and = 27°, and the line of 
sight is through the unshocked WR wind. Crel e was adjusted in 
each model to obtain a specific flux at 22 GHz. 



Figure 9. The influence of the stellar luminosities on the model 
radio spectra. Both models had n = 0.11, p = 2.0, Crcl c = 0.215, 
Cb = 2.6 X 10~ 4 , and = 50°, and the line of sight is through 
the WCR. 
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Figure 10. Model fits to radio data of WR140 at <j> = 0.837. 
In these models, which span the range of 7] considered, the low 
frequency turnover is due to free-free absorption. The data are 
represented by s o lid ci rcles, with fluxes taken from Table 2 of 
iDoughertv et all l2005l> : in this and subsequent figures we as- 
sume that the uncertainty in the a bsolute flux scale is 5 % for 
v > 15 GHz and 2% for v < 15 GHz iPerlev fc Tavlorl2003l) . The 
model parameters are noted in Table 151 and the line of sight is 
through the unshocked WR wind. Model A has r\ = 0.22, while 
model B has r) = 0.02. The low frequency turnover is shallower 
in model A, and the high frequency total emission has a stronger 
concave curvature, relative to model B. The variation of the high 
frequency emission is due to differences in the slope of the syn- 
chrotron and free-free emission in both models. The thermal flux 
from the unshocked winds and the WCR is ^ 4 mjy at 50 GHz 
in both models. 



sightlines through the unshocked O wind require p > 2.5 
for the non-thermal electron distribution to offset the syn- 
chrotron spectrum between 5-22 GHz becoming flatter as 
9 — > 90° and the free- free opacity increases (cf. Fig.|SJ. While 
such values of p are appropriate for modified shocks, the re- 
quired values of 9 are much too high comp ared to the value 
of 52° determined by the orbital solution of lDoughertv et all 
i2005l) . In fact, we find that all model fits with p > 2 are 
unsatisfactory, with either f re i, e or 9 being uncomfortably 
high. 

Since the opening angle of the WCR varies with r], 
the value of 9 required to match the low frequency syn- 
chrotron turnover varies in a similar fashion. For the models 
in Fig. 1101 9 is 5° for model A and 35° for model B (accept- 
able fits with rj between 0.02 and 0.22 can also be obtained). 
The value o f 9 from model A is inco nsistent with the 52° de- 
termined bv lDoughertv et al.l J2005T) at this phase, and while 
in model B it is still different by some considerable margin, it 
is not entirely implausible given the large uncertainty that 
exists for 9. It is anticipated that models with r\ < 0.02 
should result in fits with 9 closer to the value determined by 
IDoughertv et alJ fe005l) . but rj is not expected to be lower 
than ~ 0.01 because the mass-loss rate of the WR star then 
becomes too large. 

The slope of the high frequency synchrotron emission 
is sensitive to both p (Fig. [7J and (b (Fig. |HJ . Since there 
is almost complete degeneracy between these parameters, 
a variety of fits with different families of p and can be 
obtained (Fig. 1 1 I i . For a given rj, the value of 9 required 
to match the low frequency turnover is largely unchanged 
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Figure 11. Model fits to the radio data at <j> = 0.837, again where 
the low frequency turnover is due to free-free absorption (mod- 
els C, D, and E). Due to the sensitivity of the slope of the high 
frequency synchrotron emission to both p and £ b > it is possible to 
obtain, for a given rj, fits of similar quality with different values 
of p and B. All models had rj = 0.11 (see Table 151 . 
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Figure 12. Models of the radio data at <j> = 0.837 which demon- 
strate that it is possible to obtain fits of similar quality with a 
higher value of p when the stellar luminosities are reduced. Both 
models (B and F) had rj = 0.02 and Cb = 0.05 (see TableEJ. 



in these fits. If the stellar luminosity is reduced, fits with 
a higher value of p can be obtained, as shown in Fig. 1121 
However, we are unable to fit the data with rj = 0.02 and 
p — 2, despite making the other parameters as favorable 
as possible for this task (e.g., £b = 0.05 and luminosities 
reduced by a factor of two). 

As already noted, we are also able to obtain reasonable 
fits to the data when the Razin effect is responsible for the 
low frequency turnover. Again, fits of similar quality can be 
obtained across a wide range of rj and p, as demonstrated in 
Fig. 1131 In these models an attempt was made to match 9 as 
closely as possible to the va l ue det ermined from the orbital 
solution of lDoughertv et al.N2005l) . 

Before the relative merits of the models with either the 
Razin effect or free-free absorption being responsible for the 
low frequency turnover are considered, it is worth looking 
at the make-up of each of these models. The intrinsic syn- 
chrotron (i.e. before circumstellar absorption) and the free- 
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Figure 13. Model fits to the radio data at <fr = 0.837 where the 
low frequency turnover is due to the Razin effect (models G, H, 
and I). The model parameters are again noted in Tabic 13*1 
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Figure 14. The amount of free-free absorption affecting the in- 
trinsic synchrotron emission in model A, where the low frequency 
turnover is due to free-free absorption. Various emission compo- 
nents are shown - free- free flux (long-dashed), synchrotron flux 
(short-dashed), intrinsic synchrotron flux (dotted), and total flux 
(solid). The spectral index of the intrinsic synchrotron emission 
is —0.34 between 1.5 and 5 GHz, and steepens to —0.47 between 
15 and 22.5 GHz. This illustrates the effect of IC cooling, since, 
for p = 1.4, a = —0.2 is expected. 



free emission components to the total flux for model A where 
the low frequency turnover is due to free-free absorption are 
shown in Fig. 1141 The free-free flux is negligible in compari- 
son to the synchrotron flux, which only suffers significant 
absorption at low frequencies. In this case, the observed 
synchrotron luminosity between 1-40 GHz is roughly 90% 
of the intrinsic synchrotron luminosity. Somewhat paradox- 
ically, the intrinsic synchrotron emission from models where 
the low frequency turnover is caused by the Razin effect can 
experience more severe free-free absorption, as seen from a 
comparison of Figs. I14| and 1151 This is because the free-free 
absorption is smallest when the sight lines are closely par- 
allel to the WR shock (as is the case for model A), and is 
significantly greater when the sight lines are deeper within 
the WCR (e.g., compare the 10 GHz flux in Fig. when 
6 = 10 or 20° with the flux when 6 = 40° - see also the dot- 
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Figure 15. As Fig. lMl but for model J, where the Razin effect is 
responsible for the low frequency turnover. 



Table 3. Parameters for model fits to the <j> = 0.837 radio data. 
In models A-F the low frequency turnover is due to free-free ab- 
sorption, while in models G-J it is due to the Razin effect. The 
values of for which the sight lines are parallel to the asymp- 
totic opening angle of the WR shock, contact discontinuity, and 
O shock are noted in Table |21 



Model 


V 


V 


Crel,e 




Cb 




e 


A 


0.22 


1.4 


1.38 x 10" 


3 


0.05 




5° 


B 


0.02 


1.4 


5.36 x 10" 


3 


0.05 




35° 


C 


0.11 


1.53 


2.26 x 1(T 


4 


0.5 




13° 


D 


0.11 


1.4 


2.03 x 10- 


-3 


0.05 




13° 


E 


0.11 


1.1 


2.79 x 10- 


2 


5.0 x 10- 


-3 


14° 


F 


0.02 


1.6 


1.81 x 10" 


-3 


0.05 




40° 


G 


0.11 


2.0 


0.22 




2.6 x 10- 


4 


50° 


H 


0.11 


1.6 


0.15 




4.0 x 10- 


4 


44° 


1 


0.0353 


1.4 


0.14 




1.0 x 10- 


-3 


58° 


,) 


0.22 


1.5 


0.14 




4.0 x 10- 


4 


36° 



ted curve in Fig. 12b in lPittard et alj|200fj) . The observed 
synchrotron luminosity in Fig. 1151 is only 40 per cent of the 
intrinsic synchrotron luminosity, and in such cases the IC 
luminosity predicted using Eq. [3] will be substantially un- 
derestimated if the absorption is not accounted for. 

4-3.3 Synthetic radio images 

The spatial distributions of radio emission from models A 
and B are shown in the left and right panels of Fig. 1161 
In all panels the intensity of emission from the WCR ex- 
ceeds that from the stars, though the relative amount of 
free- free emission increases strongly with frequency. Model A 
is viewed at low inclination, while model B is viewed at a 
somewhat higher inclination (Table |3J. The absorption of 
emission from the far side of the WCR by the O wind is 
negligible in model A, but significant in model B, as evident 
by the dearth of emission near the O star. As expected, this 
hole is much larger at 1.7 GHz than at 8.4 GHz. The O star 
is not located at the centroid of the hole, since the O wind is 
restricted on one side by the WCR. Those parts of the WCR 
with optically thin lines of sight through the outer envelope 
of the "modified" O wind are predominantly observed. In 
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Figure 16. Intensity distributions at 1.7 GHz (top) and 8.4 GHz (bottom) from model A (left) and model B (right) of WR140. The 
contour levels are spaced apart by 0.5 dex, are chosen to highlight the large-scale structure, and are identical in both model's at a specific 
frequency, though change with frequency. The positions of the stars are indicated by crosses (WR star to the upper-right of the O star). 
The effect of different wind momentum ratios and line of sight angles are clearly visible. 



addition, the intensity of the far side of the WCR is highest 
nearer the apex. The reduced opening angle of the WCR in 
model B relative to model A is clearly visible. IC cooling 
creates a "V"-shaped intensity distribution in model A, but 
this is not so clearly seen in model B, largely due to the 
difference in viewing angle. 

Though not shown here, another important finding is 
that the relative position of the stars and the observed peak 
of the emission from the WCR is not necessarily related to 
r), due to the occultation of the WCR by the O star stellar 
wind. In fact, as 9 increases and D decreases the absorption 
can become so large that the position of the WCR apex 
appears closer to the WR star than to the O star, despite 
the WR star having the stronger wind. Finally, we note that 
no offset between the center of emission at 1.7 GHz and 
8.4 GH z is seen, contrary to the report bv lDoughertv et alJ 

teoorJ) . 

To more closely compare the model intensity dis- 
tributions at 8. 4 GHz to the VLB A images shown in 
iDoughertv et all l)200fJl . we have used the AIPS subrou- 
tine UVCON to generate visibilities appropriate for a VLBA 
"observation" of our models. System noise estimates were 



accounted for by including the performance characteristics 
of the VLBA telescopes, e.g., antenna efficiencies, system 
temperatures etc. The resulting visibilities are then imaged 
and deconvolved using standard techniques. The resulting 
simulated observations are shown in Fig. 1171 With perfect 
data i.e. no noise, it is possible to distinguish the two dif- 
ferent models (Fig. I17b .'l . However, once appropriate noise 
has been incorporated, the situation is quite different as can 
be seen in Fig. 117b and c. Adopting the s ame imaging pa- 
ramet ers (imaging weights etc.) as used in IDoughertv et "all 
1)20051) . the differences between the models are not apparent. 
However, applying a more extreme taper to strongly empha- 
sise the lower spatial frequency content of the visibility data 
(Fig. 117b ) the different models can just be discerned by an 
apparent difference in the opening angle. These synthetic 
images show that any curvature in the WCR, from either 
free-free absorption or from the wind-momentum ratio, is at 
best difficult, and may be impossible, to detect with VLBI 
observations. It is clearly dependent on the spatial frequency 
content and the attainable signal-to-noise ratio of the data. 
The high-brightness emission that is detected by VLBI ar- 
rays arises predominantly in the vicinity of the shock apex, 
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where the opening angle of the shocks have yet to reach 
their asymptotic values. Even when the lower surface brigh- 
ness emission from electrons further down the post-shock 
flow can be detected, e.g., by arrays with shorter baselines 
to recover lower spatial frequency visibilities, the beam of 
the array smooths out the features of the WCR. This leads 
us to call into question the reliability of values of rj derived 
from these types of observation. 

4-3.4 Discussion of the radio models 

In attempting to identify the most suitable models, it is 
noted that the fits in Fig. llHI reauire a low B-field within the 
WCR (£b ~ 3 x 10 -4 — 10 -3 ), which is far from equipartition 
with the energy density of non-thermal particles (electrons 
plus ions). Assuming a simple dipole field, this implies that 
the B-field at the surface of each star is very low, with little 
or no amplification in the WCR. As the synchrotron flux is 
proportional to frel.e Cb^ 4 ' l ar S e values of frei.e are required 
to normalize the model spectra to the data. The models 
(G-J) shown in Figs. IT31 and IT51 imply that 10-20 % of the 
available wind kinetic power processed through the shocks 
bounding the WCR (which is a small fraction of the total 
kinetic power of the winds in the system) is transferred to 
the non-thermal electrons. This is unacceptably high. Hence, 
models where the low frequency turnover is due to free-free 
absorption are preferred. The high energy non-thermal emis- 
sion from models G-J is also less consistent with current 
EGRET and INTEGRAL observations, and in model G the 
non-thermal X-ray flux exceeds the observed thermal flux. 
This is discussed further in Sec. 14.41 

Amongst models A-F, model E is notable for an ex- 
tremely low value of p. The resulting high energy non- 
thermal spectrum is too flat if WR 140 is indeed associated 
with the EGRET source 3EG ,12022+4317 (see Sec. K£2\ . 
The magnetic energy density in model C is towards the 
high end of expectations, with B = 3.3 G at the apex of 
the WCR. A magnetic field this high is likely to prevent 
particle acceleration as the shock velocity does not exceed 
the phase velo city of whistler waves propagating normal to 
the shock (see lEichler fc Usovlll993l for the necessary con- 
ditions) . At this stage we also prefer to use th e stellar lumi- 
nosities determined in lDoughertv et alj (l2005f) . rather than 
reduced values such as used in model F, though further in- 
vestigation of this matter is needed. On this basis, models 
A, B and D are preferred. Of these three models, the view- 
ing angle of model B i s closer to the value determined by 
iDouehertv et all <l2005l) . and hence is the preferred model 
overall. 

We now compare our estimate of rj from model B to 
previous estimates. First, rj — 0.02 i s con sistent with the 
value derived from iMarchenko et al.l d2003T) . from an anal- 
ysis of the CHI and H el line profiles around periastron. 
IDouehertv et all <l2005t) estimate a value an order of mag- 
nitude higher (0.22), but this is an indirect estimate based 
on an assumed mass-loss rate for the O supergiant. They 
also argue that rj — 0.22 is consistent with the open- 
ing angle of the WCR as estimated from the VLBA im- 
ages. However, as argued in the previous section, we should 
not expect the opening angle of the WCR to be con- 
strained by the VLBA observation s at this orbital phase. 
IVarricatt. Williams fc AshokH2004l) . from Hel absorption in 
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Figure 17. 8.4 GHz synthetic VLBA observations of the model 
intensity distributions in Fig. 1161 The data have been scaled 
to ma tch the peak flux in the observations in fPoughcrt v et all 
i2005l) at phase 0.841 - close to 1.2 mjy beam -1 . The top panel 
(a) shows the resulting image for model A (thin) and model B 
(thick) when there is no noise i.e. "perfect" data. The middle 
panel (b) used the same ima ging and deconvolution parameters 
as in [Dougherty et all 120051) , and shows no discernable differ- 
ence in the images of the two models. In the lower panel (c), the 
lowest spatial frequency content of the data is more strongly em- 
phasized by severely tapering the data, and the two models may 
just be distinguished. Contours in (a) start at 30 fijy beam -1 
and increase by powers of 1.5. For (b) and (c), the first contour 
is at 33 jijy beam -1 . Each minor tick on the axes represents 
1 mas. The scale and distribution of the emission in panel (b) is 
consistent w ith the observations at both phase 0.830 and 0.841 
presented in lDoughertv et al.l J2005fl . 
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the WR wind, determine a lower limit to the half-opening 
angle of 42°, which is consistent with the value for the WR 
shock in model B. 



4.4 The high energy emission 

4-4-1 Models of the high energy emission 

The parameters obtained from the models of the radio 
data in Sec. 14.31 allow prediction of the IC and relativistic 
bremsstrahlung flux, and by specifying the ratio of the en- 
ergy density of electrons to ions it is possible to predict the 
flux from the decay of neutral pions. Such calculations can 
be used to further constrain the spectral index of the non- 
thermal paricles and thus the nature of their acceleration, 
as demonstrated in this section. 

The predicted non-thermal emission from model B is 
shown in Fig. 1181 and a summary of the emission compo- 
nents noted in Tabled The INTEGRAL flux is insensitive to 
reasonable assumptions about the value of 7 max , though the 
EGRET and GLAST fluxes are more so. However, the VER- 
ITAS flux is highly sensitive to both 7 ma x and Crei,i/Crei,e, 
and the predictions presented here are meant only as a guide. 
At GeV energies and above we may underestimate slightly 
the flux through our neglect of the emission from secondary 
particles. 

The total luminosity from non-thermal processes is 
~ 0.5% of the kinetic power in the wind-wind collision 
(~ 10 36 ergs -1 ), the latter being ~ 1% of the total kinetic 
power in the winds. The IC emission dominates for photon 
energies less than 50 GeV, while the emission from pion de- 
cay reaches energies up to 15 TeV (the maximum photon 
energy from each process is limited by the assumed value 
of 7max = 10 5 for the electrons and ions). The relativistic 
bremsstrahlung emission is a minor contributor to the to- 
tal non-thermal flux, and does not dominate at any energy. 
The absorption of high energy photons by two-photon pair 
production is significant at GeV and TeV energies, in good 
agreement with the calculation in Fig. [5] where a single line 
of sight from the stagnation point of the WCR to an observer 
was considered. 

The IC flux is directly proportional to the energy den- 
sity of non-thermal electrons, U le \,e, and thus to Crei.e, which 
in the calculation for Fig. llBl was set to 5.36 x 10~ 3 as deter- 
mined from the fits of the radio data (Fig. 1101 . The gradual 
steepening of the IC spectrum between 2-100 keV (corre- 
sponding to 14 < 7 < 100) indicates the effect of IC cooling 
on the non-thermal electron energy spectrum. The lack of a 
clearly defined spectral break reflects the situation for the 
sy nchrotron spectrum at radio frequencies (see Figs. 4 and 6 
in lPittard et alJl2006t) . Above 1 MeV, the spectral index is 
equal to — p/2 = —0.7 (corresponding to a photon spectral 
index, r = -a + 1 = (p + 2)/2 = 1.7). 

The relativistic bremsstrahlung spectrum is fairly flat 
below 1 MeV, and falls off sharply at higher energies. Like 
the IC emission, its flux is directly proportional to the spec- 
ified value of Crei.e- In contrast, the emission from neutral 
pion decay peaks at 0.5m„-c 2 = 67.5 MeV, and, for pions 
produced from collisions involving non-thermal protons, has 
a maximum energy of -E p , m ax/6, where the maximum en- 
ergy of non-thermal protons, -E p ,max = 7maxm p c 2 . The flux 
from pion-decay is directly proportional to mU Te i t iV, where 



m is the number density of thermal ions, f7 re i,i = Crei.iUth 
is the energy density of non-thermal ions, and V is the 
volume occupied by the WCR. Since U t h oc m oc D~ 2 , 
Urel,i oc £ r ci,i-D~ 2 . In the adiabatic limit, the volume of the 
WCR scales as D 3 , so the intrinsic flux from neutral pion 
decay should scale as Crei,i.D - . As we do not know how (or 
if) Crci.i varies with D, the overall dependence on D is un- 
known. However, since the optical depth from two-photon 
pair production increases as D^ 1 , the observed flux above 
100 GeV is likely to decline as D is reduced. An added com- 
plication is that 7 m ax for the non-thermal ions will increase 
with decreasing D if the stars move sufficiently close that 
the geometry of the magnetic field in the unshocked stellar 
winds changes from toriodal to radial. In the radial limit, we 
expect that 7 ma x oc 1/D, a nd thus 7 ma x m ay reach values 
of 10 6 in close systems Ccf. lBednarekll2005l) . Table 0] sum- 
marizes the emission from the model components in various 
instrumental bands. 

Fig. 1181 also compares the non-thermal X-ray and 7- 
ray emission calculated from our model with the observed 
ASCA data at cj> = 0.837 - dataset 27022010, observed on 
1999-10-22. The O star is in front of the WR star at this 
orbital phase, and although the IC emission will be slightly 
higher when a full treatment of the anisotropy of the scat- 
tering is accounted for, t he increase should be less tha n a 
factor of 3 (see Fig. 10 in lReimer. Pohl fc Reimeill2006l) . It 
is reassuring that even with this factored in, the model IC 
fluxes are significantly lower than the observed X-ray emis- 
sion, consistent with the lack of detection of an underlying 
power-law component. The inset of Fie. ll8l shows the intrin- 
sic and attenuated thermal X-ray emission from model B. 
An important point is that the attenuated emission is in 
close agreement with the ASCA data. Unfortunately, this is 
also the case for our model A, and so the X-ray absorption 
at phase 0.837 cannot be used to constrain 77. 

Models of the radio data alone are ill-constrained (see 
Sec. 14. 3L and the properties of their high energy non- 
thermal emission provides an important tool for discriminat- 
ing between the models. In Fig. 1191 . the high-energy spectra 
of models A and B, are shown. The flux from model B is 
slightly higher than from model A because f re i, e and £ re i,i 
are slightly higher, and the IC emission has a slightly steeper 
spectral slope in model B due to stronger IC cooling of the 
non-thermal electron energy spectrum resulting from the 
smaller distance between the apex of the WCR and the O 
star. 

Fig. 1201 shows the high energy non-thermal emission 
from models C, D, and E. The degeneracy of the emission at 
radio wavelengths is broken at high energies where clear dif- 
ferences in the models occur. In particular, the decrease in p 
from model C (p = 1.53) to model E (p = 1.1) is clearly man- 
ifest in the spectral slope, while the increase in £rei,e from 
model C to E is apparent in the relative normalization of 
the spectra. Since model E has £ re i ?e = 2.8 x 10~ 2 which pre- 
cludes setting f re i,i = 100Crel,e, we instead set £ re i,i = 30£roi,o 
for this model. 

The high energy non-thermal emission from models B 
and F are compared in Fig. 1211 where the main difference 
is in the assumed stellar luminosities. The spectral slope of 
the IC emission at energies greater than 0.1 MeV reflects the 
different values of p in these models. This emission arises pre- 
dominantly from those electrons near the apex of the WCR 
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Figure 18. The radio, and non-thermal UV, X-ray and 7-ray emission calculated from model B, together with the observed radio and 
X-ray flux. The radio model shown indicates the thermal free- free flux (displayed below the data points), the intrinsic synchrotron flux 
before free- free absorption (displayed above the data points), and the total observed emission (cf. Fig. I1UI . The model IC (long dash), 
rclativistic bremsstrahlung (short dash), and pion decay (dotted) emission components are shown, along with the total emission (solid). 
Photoelectric absorption is significant at soft X-ray energies up to ~ 2 keV. To calculate the IC and relativistic bremsstrahlung emission, 
Crel e = 5-36 X 10 -3 , as derived from the fit to the radio data (Fig. ITTJI was assumed. Crcl , i = 100 X £ re i e was assumed to calculate 
the neutral pion decay emission. The observed ASCA X-ray spectrum from <j> = 0.837 is also shown. The intrinsic (dashed line) and 
attenuated (dotted line) thermal X-ray emission from model B is displayed in the inset, where we have assumed unmodified collisional 
shocks, rapid thermalization of the post-shock electrons, and collisional ionization equilibrium in our model. The ISM absorbing column 
was fixed at TVjj = 5.4 X 10 cm" " 2 IZhekov fc Skinnerll200ri). and the ab sorption coefficients for the ISM and the WR and O star winds 
were calculated using Cloudy v94.00 iVan Hoof. Martin fc Ferlandll2000l see |http: / / www.nublado.org) . 



Table 4. The high energy non-thermal emission from model B of WR 140 at orbital phase 0.837, with n = 0.02 and f re i e = 5.36 X 10 — 3 . 
Crel i = 100 X Crel e i s assumed for the calculation of the emission from neutral pion decay. The IC emission dominates the total emission 
in each energy band noted in this table. 



Mechanism 


Total emission 


EGRET 


EGRET 


INTEGRAL IBIS 


GLAST 






(lOOMeV-lOGeV) 


(lOOMeV-lOGeV) 


(15keV-10MeV) 


(20MeV-300GeV) 




(ergs" 1 ) 


(ergs^ 1 ) 


(photons s _1 cm - 2 ) 


(photons s— 1 cm— 2 ) 


(photons s— 1 cm -2 ) 


Inverse Compton 


3.5 x 10 33 


2.1 X 10 33 


5.4 X 10~ 9 


2.7 X KT 6 


1.5 x 10- 8 


Rel. bremsstrahlung 


1.8 x 10 3() 


3.2 x 10 29 


2.5 x 10- 12 


2.0 X 10" 9 


2.1 X 10~ n 


ir° decay 


6.8 x 10 31 


9.8 x 10 30 


1.3 x 10- 11 


3.9 x 10- 14 


1.8 x 10" 11 



with initial post-shock Lorentz factors greater than ~ 100. 
All such electrons rapidly cool to Lorentz factors of approxi- 
mately 12 (model B) and 25 (model F) before advecting out 
of the system. Although the absorption of 7-rays by pair- 
production is lessened in model F due to the reduced stellar 
luminosities assumed, the commensurate reduction in £ re i >e 
(and thus Crei.i) and the increase in p lead to an overall de- 
crease in the VERITAS flux (Table 0. 

Fig. 1221 shows a comparison between the high energy 



non-thermal emission from models where free-free absorp- 
tion (model B) or the Razin effect (models G, H, I and J) is 
responsible for the low frequency radio turnover. The large 
values of C re i, e in models G-J cause the IC flux to be at 
least two orders of magnitude higher relative to model B. 
Future observations with AGILE and GLAST should have 
the sensitivity to detect such fluxes, thus allowing a pro- 
cess of discrimination between our models. In fact, we can 
already exclude model G, since its predicted non-thermal 
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Figure 19. The high energy non-thermal emission calculated 
from models A and B with p = 1.4. 
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Figure 21. The high energy non-thermal emission calculated 
from models B and F, with 77 = 0.02. 





Figure 20. The high energy non-thermal emission calculated 
from models C, D and E, with r\ = 0.11. 



Figure 22. The high energy non-thermal emission calculated 
from models B, and G-J. The low frequency radio turnover in 
model B is due to free-free absorption, while in models G-J it is 
due to the Razin effect. The ASCA data are also shown. 



emission exceeds the observed thermal emission at X-ray 
energies. While models H-J cannot be excluded in the same 
manner, there are other reasons why they are less desirable 
than model B (see Sec. 14. 3, 21 . Table|S]notes the non-thermal 
luminosities and fluxes in various energy bands from each 
model. 

4-4-2 Comparison to observations 

The total photon flux in the EGRET band from 
models A-F (where free-free absorption is responsi- 
ble for the low frequency radio turnover) varies from 
4.8 x ICT 10 photons s^cm -2 (model C) to 4.2 x 
10~ 8 photons s -1 cm~ 2 (model E). The highest prediction is 
still below the flux detected from 3EG J2022+4317 during 
the December 1992 observation at (f> = 0.97 (2.47 ± 0.52 x 
10 -7 photons s _1 cm -2 ). While the different orbital phases 
of the radio and EGRET data prevent a direct comparison, 
the IC emission should increase between phase 0.837 and 
0.97 as the stellar separation decreases by about a factor of 
three. Thus, the models with the highest predicted fluxes 
are broadly consistent with the data. In contrast, models G, 
H and J have predicted photon fluxes in the EGRET band 



which exceed that observed from 3EG J2022+4317. This is 
an additional reason why these models are not as desirable 
as those where free-free absorption is responsible for the low 
frequency radio turnover. 

In model B, the photon spectral index of the total non- 
thermal emission, V, is -2.9 at 1 keV (due to photoelectric 
absorption), 1.7 at 1 MeV and 100 MeV, 1.8 at 1 GeV, and 
2.5 at 10 GeV, and is in rough agreement with the observed 
EGRET value of V = 2.31T0.19 between 100 MeV- 10 GeV. 
Model E has the flattest non-thermal spectrum, with F = 
-1.0 at 1 keV, 1.55 at 1 MeV and 100 MeV, 1.6 at 1 GeV, 
and 2.1 at 10 GeV, and matches the observed spectral slope 
less well. 

4-4-3 Comparison to previous work 

The non-thermal 7-ray e mission from WR 140 has p revi- 
ously been estimated by iBenaglia fc Romerol feOOSt and 
iReimer. Pohl fc Reimed (12006^1 . Since Benaglia & Romero's 
calculations are based on a stellar separation of 10 AU (cor- 
responding to (f> = 0.95), neither our estimates, their esti- 
mates, or the EGRET detections, can be directly compared. 
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Table 5. The high energy non-thermal emission calculated from the models in Table 131 



Model 


Total emission 


EGRET 


EGRET 


INTEGRAL IBIS 


GLAST 


VERITAS 






(lOOMeV-lOGeV) 


(lOOMeV-lOGeV) 


(15keV-10MeV) 


(20MeV-300GeV) 


(200GeV-50TeV) 




(ergs -1 ) 


(ergs" 1 ) 


(photons s" 1 cm -2 ) 


(photons s" 1 cm -2 ) 


(photons s~ 1 cur 2 ) 


(photons s" 1 cm" 2 ) 


A 


2.4 X 10 33 


1.5 x 10 33 


3.4 x 10 -9 


1.0 X 10~ 6 


1.0 X 10~ 8 


2.7 X 10~ 14 


B 


3.6 X 10 33 


2.1 x 10 33 


5.4 x KT 9 


2.7 X 10" 6 


1.5 X 10~ 8 


4.5 X 10" 14 


C 


3.2 x 10 32 


1.8 x 10 32 


4.8 x 10 -10 


3.0 x 10~ 7 


1.4 x 10" 9 


3.5 x 10" 15 


D 


2.8 x 10 33 


1.6 x 10 33 


4.0 x 10" 9 


1.5 x 10" 6 


1.1 x 10~ 8 


2.8 x 10" 14 


E 


3.6 x 10 34 


2.0 x 10 34 


4.2 x KT 8 


4.2 x KT 6 


9.4 x 10" 8 


2.8 x 10" 13 


F 


1.3 x 10 33 


7.0 x 10 32 


2.0 x KT 9 


2.4 x KT 6 


6.3 x 10~ 9 


1.5 x 10" 14 


G 


2.4 x 10 35 


7.6 x 10 34 


2.8 x 10~ 7 


1.4 x 10" 3 


1.1 x 10~ 6 


8.0 x 10" 15 


11 


1.9 x 10 35 


1.1 x 10 35 


3.1 x 10~ 7 


2.7 x 10~ 4 


9.5 x 10" 7 


9.9 x 10" 15 


1 


1.1 x 10 35 


6.8 x 10 34 


1.7 x 10" 7 


8.1 x 10" 5 


4.7 x 10" 7 


6.5 x 10" 15 


.1 


2.4 x 10 35 


1.5 x 10 35 


4.0 x 10" 7 


1.7 x 10" 4 


1.1 x 10" 6 


1.1 x 10" 14 



Their estimate for the IC luminosity in the EGRET band of 
2.1 x 10 34 erg s -1 is typically higher than the values from our 
models A-F, but lower than those from models G-J (cf. Ta- 
ble EJ. 

Th ere are substantial dif f erence s between this work and 
that of iBenaglia fe Romerol ll200otl : their estimate is ob- 
tained through the use of Eq. [3J in which the ratio of the 
IC to synchrotron emission depends strongly on the B-field, 
and is based on only a single position in the WCR. They es- 
timate that B = 200 mG (although noting its value is quite 
uncertain). In contrast, here the IC emission is determined 
from models of the radio emission. The magnetic field is a 
fitted (though ill-constrained) parameter of the models, with 
a maximum at the stagnation point 3 , and rapidly declines 
in the downstream flow. Also, in the approach presented 
here the B-field directly affects the synchrotron emission, 
but also indirectly affects the calculations of the high energy 
emission which depend on the population and spatial distri- 
bution of non-thermal par ticles i.e. £ re i, e and Crci.i- Finally, 
IBenaglia fe Romerol i2005t) use the observed synchrotron lu- 
minosity, rather than the intrinsic luminosity, leading to an 
underestimate of the IC luminosity. 

Our estimates for the relativistic bremsstrahlung lumi- 
nosity ar e broadly consistent with th e limit of < 10 32 erg s _1 
noted in IBenaglia fe Romerol teOO&tl . but the pion decay lu- 
minosities are wildly discrepant with our lowest estimate 
being 10 orders of magnitude high e r. Det ails of their calcu- 
lation are noted in IBenaglia et afl i200lf) . where a "cosmic 
ray enhancement factor" of 10 is assumed. This factor is the 
ratio between the cosmic ray energy density at the source (in 
this case the WCR) and that observed at the Earth (which 
is roughly equal to the average cosmic ray energy density 
within the Galaxy of 1 eVcm~ 3 = 1.6 x 10~ 12 erg cm" 3 ). 
However, at the head of the WCR, the energy density of rela- 
tivistic ions in our model B, for example, is ~ 0.64 erg cm -3 . 
This corresponds to a cosmic ray enhancement factor of 
~ 10 11 , and accou nts for the large underestim ate of the pion 
decay emission by IBenaglia fc Romerol j200ot) . 



3 Since the magnetic field is not well constrained by the radio 
fits, there is a large variation in its strength between models. For 
instance, at the apex of the WCR, B = 74 mG in model G, rising 
to 3.3 G in model C, and is 1.2 G in our preferred model (B). 



iReimer. Pohl fc Reimerl (120061) have also calculated the 
expected high energy non-thermal flux from WR 140. They 
conclude that the change in the IC flux with viewing an- 
gle due to anisotropic scattering is likely to be obscured by 
large variations in the energy density of the stellar radia- 
tion fields resulting from the high orbital eccentricity. While 
anisotropic IC emission and Klein-Nishina effects are clearly 
processes which should be included in future models, our 
work shows that other uncertainties, such as the value of p, 
currently have a greater influence on the IC emission. 

4-4-4 Predictions for forthcoming instruments 

Our predictions for the flux in the energy band of the IBIS 
instrument onboard INTEGRAL and the GLAST large area 
telescope (LAT) are also noted in Table [S] For models A-F 
the predicted fluxes are ~ 10 3 x lower in the INTEGRAL 
band, a nd ~ 10 2 x lower in the G LAST band than those 
made bv lBenaglia fc Romerol <200ot) . Using the INTEGRAL 
Observing Time Estimator 4 , we find that model A cannot be 
detected at the 5a level in the 15 keV— 10 MeV range using 
the IBIS imager, since an exposure time of 300 yrs is re- 
quired! Model E, the brightest of models A-F in this band, 
still requires an exposure of 180 yrs. With several bright 
sources (including Cyg X-l) nearby, the required exposure 
time may be even longer. On the other hand, for model J 
an exposure time of only 2.8 Msec is needed. At the time 
of writing, WR 1 40 is not detected in a 2 Msec exposure 
JPe Beckeill2005lh and we are close to ruling out model J 
(and H) on this basis. Line emission at 7-ray energies from 
the de-excitation of nuclear isomers, such as 12 C* and 16 O*, 
is m uch too faint to be detecte d with current instrumenta- 
tion iBenaglia fc RomerolEool . 

The GLAST 5a sensitivity at E > 100 MeV for sources 
at high galactic latitudes with a F = 2 spectrum after a 
2 yr all-sky survey is 1.6 x 10~ 9 photons s" 1 cm -2 . There- 
fore, our predictions suggest that WR 140 will be detectable, 
this being important to differentiate between the models in 
Table 

At energies of approximately 50 GeV and higher, the 
flux from the models arises exclusively from the decay of 

4 http:/ /integral. esac.esa.int/isoc/operations/html/OTE. html 



Models of the Colliding Winds Binary WR140 23 



neutral pions, with a high energy cut off of « 15 TeV, 
strongly dependent on 7 max of the ions. The VERITAS- 
4 array is sensitive to 7-rays in the 200 GeV to 50 TeV 
energy band, but unfortunately the high absorption from 
two-photon pair production severely reduces the photon 
flux above ~ 10 GeV. The predicted flux from model B 
is not sufficient for a detection in a 50 hr observation 
(cf. iFeean. Hall fc Vassilievl l2003h . and we conclude that 
WR 140 is probably too faint to be detected at <j> = 0.837, 
though it might be brighter at phases closer to periastron. 
An observation by GLAST, for example, will better con- 
strain the potential models, and allow more robust predic- 
tions of the flux at ~ TeV energies. 



5 DISCUSSION 

5.1 Mass-loss rates of early- type stars 

It is important that accurate measurements of mass-loss 
rates from early-type stars are obtained, as mass-loss is 
known to significantly affect their evolution, and the injec- 
tion of mass, momentum and energy influences not only the 
gas-phase conditions in the immediate environment of the 
clusters within which the massive stars form, but also the 
phase structure and energetics of the ISM on galactic scales, 
and the thermodynamics and enrichment of the intergalactic 
medium on Mpc scales and above. 

Unfortunately, observational determinations of mass- 
loss rates are often highly uncertain, in part because many 
of the techniques are sensitive to the degree of structure 
in the winds. Overestimates by factors of at least 3 or 
more are likely when the winds are clumpy, but for some 
types of O stars it could be as high as factors of 20 or 
more iFullerton. Massa' fc Prinial lioOrj . The X-ray emis- 
sion from CWB systems provides an opportunity to mea- 
sure mass-loss rates using a method which is not expected 
to be very sensitive to clumping (see Sec. 14.2.31 . In fact, 
mass-loss rates could be underestimated using this tech- 
nique: at least so me winds may be both denser and faster 
at high latitudes JPwarkadas fc Owockill2002l; ISmithll2002l 
IChesneau et al.l l2005rl . and if the stellar rotation axes of 
the stars are closely perpendicular to the orbital plane, the 
WCR then probes those parts of the winds which are less 
dense and slower. In addition, a significant amount of energy 
may be siphoned into particle acceleration. In this work, 
the mass-loss rate for the O supergiant star in WR 140 is 
towards the lower end of previous estimates, but is consis- 
tent with expected biases. A reduction in mass-loss rates, 
and/or mass-loss through a predominantly polar flow, can 
also explain the near-symmetry of observed X-ray lines from 
early - type stars dOwoc ki fc Cohenl2006t iMullan fc Waldronl 
l2006t ICohen et alJ^OodT and both reduced mass-loss and 
polar fl ow help in the theoreti cal formation of gamma-ray 
bursts dWooslev fc Hegerll2006ft . 

5.2 The nature of shocks and particle acceleration 
in wide CWB systems 

In this work we have presented evidence for a fairly hard 
spectrum (i.e. p < 2) of non-thermal electrons in WR 140, 
which is further supported by earlier fits to the radio data 



of WR 147 (this was mistakenly interpret ed as evidence for 
shock modification in lPittard et al.ll2006l) . As discussed in 
Sec. 12.1.31 there are many potential mechanisms for pro- 
ducing p < 2. One possibility is particle re-acceleration, 
but while the wind-embedded shocks believed to occur in 
stellar winds due to the line de-shadowing instability o f 
the radiative driving ('e.g.. lOwocki. Castor fc Rvbickilll988f) 
may provide such an opportunity, and thus harden the 
non-thermal particle spectrum prior to the global shocks 
bounding the WCR, further evidence is needed for this 
explanation to be more persuasive, not least because re- 
acceleration has yet to be demonstrated in theoretical cal- 
culations of wind-embedded shocks. One problem is that 
in current one- dimensional simulations of wind-embedded 
shocks, the accelerated particles are trapped between for- 
ward and reverse shocks. This has posed a problem for 
the in terpretation of non-thermal radio em ission from single 
stars llVan Loo. Runacres fc Blommdl2006l) . and only calcu- 
lations using a phenomenological sho ck model are currentl y 
successful in fitting the data (e.g., IChen fc White! Il994h . 
Two- and three-di mensional models of instab ility-generated 
structure (see, e.g.. iDessart fc OwockilfeoOS!) may allow ac- 
celerated particles to escape one shock and be re-accelerated 
at others, but this has still to be demonstrated. 

We can speculate about additional effects which occur 
in CWB systems with highly eccentric orbits. For example, 
if re-acceleration is important, one might expect p to de- 
pend on the orbital phase, as the number of wind-embedded 
shocks that accelerate a non-thermal particle prior to the 
WCR may change significantly as the stellar separation al- 
ters. Similarly, if the obliquity of the shocks changes with 
stellar separation, the injection efficiency and rate of parti- 
cle acceleration may also alter. In addition, it is likely that 
the maximum energy obtained by electrons and ions varies 
with stellar separation. Finally, the assumption of isotropic 
synchrotron emission will need to be reconsidered if the mag- 
netic turbulence is low enough that the mean B-field has a 
preferred direction. In the context of WR 140, these are all 
potential explanations for the observed asymmetry of the 
radio lightcurve. 

Is the efficiency of particle acceleration in WR 140 high 
enough for shock modification to occur? In Tycho's SNR, 
the concave curvature of synchrotron radio emission is clear, 
though sti l l indi rect, evidence for shock modification (see 
IVolk et al.ll2002l) . Unfortunately, in WR140 the observed 
synchrotron emission extends over a much narrower fre- 
quency range and a concave curvature cannot easily be dis- 
tinguished. And while the lower than expected thermal X- 
ray temperature is consistent with shock modification, a 
more probable cause is T c < 2]. Likewise, while shock mod- 
ification may help to explain the behaviour of the X-ray 
lightcurve which does n ot follow the expected 1/D variatio n 
for an adiabatic WCR JPollock. Corcoran fc StevenslEool) . 
there is a realistic chance that radiative cooling is the cause. 
Finally, the B-field in WR 140 is sufficiently high to cause 
significant energy to be transferred from the non-thermal 
particles into magnetic turbulence via Alfven wave heat- 
ing, thereby lowering the acceleration efficiency, although 
it remains possible that shock modification occurs predom- 
inantly through particle loss (due to the curvature of the 
shocks) rather than high acceleration efficiencies. To summa- 
rize, there is currently no strong evidence for shock modifi- 
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cation in WR 140, and on balance we conclude that strongly 
modified shocks are unlikely to occur in WR 140. 



5.3 Non-thermal X-ray, TeV and muon neutrino 
emission from CWBs 

While there is no observational evidence for non-thermal X- 
ray emission in WR 140, there are some systems that appear 
to show a power-law tail. Perhaps the strongest candidate 
is rj Carinae, where X-rays with a photon index A « 2.4 
were detected at ener gies of up to 150 keV during a long 
BeppoSAX exposure iViotti et alJl2004T) . Determining the 
spectal slope at higher energies should allow discrimination 
of the competing models proposed to explain this emission. 

The lack of observed non-thermal X-ray emission from 
CWB systems which ar e known non -thermal radio emitters 
<De Becker et alJl2004al.l2005ll20o3l contrasts with the pos- 
sible detection of non-thermal X-ray emission from the 3.4 d 
orbital period system HD 1 59176, which is not know n to be a 
non-thermal radio emitter llDe" Becker et al l2004bl) . Though 
the nature of the wind-wind collision in such short period 
systems is uncertain, if particle acceleration occurs, perhaps 
mainly through reconnection, IC losses will be catastrophic, 
and it is not surprising that non-thermal X-ray emission is 
possibly detected despite no clear signature of s ynchrotron 
emiss ion (due to the obvious high absorption - iDe Beckerl 
l2005h . 

In the models of WR 140 presented in this paper, ab- 
sorption of high energy photons by pair-production cre- 
ates a very steep photon spectrum between 10 GeV and 
1 TeV. Consequently, the photon flux above 1 TeV is many 
times lower than the flux predicted in the VERITAS band- 
pass, which is dominated by photons with energies less 
than 1 TeV (see Table |SJ . For example, in model B the 
flux above 1 TeV is 1.2 x 10~ 14 photons s" 1 cm~ 2 . Never- 
theless, the flux from CWB systems compares favourably 
to some other mechanisms proposed for the genera- 
tion of TeV emission, such as the interaction of cos- 
mic rays in the innermost parts of the winds of early- 
type stars llTorres. Domingo-Santamaria fc Romerol l2004h 
or within the cluster wind of a dense stellar cluster 
iDomingo-Santamaria fc TorresH2006T) . Therefore, it is pos- 
sible that several CWBs or multiple colliding winds from 
a cluster of stars contribute part or all of the flux from 
the unidentified TeV sou rce detected in t he Cy gnus re- 
gion, TeV J2032+4130 jAharonian et alJ l2005al) . which 
has an integral flux F 1 {E 1 > 1 TeV) = 6.9(±1.8) x 
10~ 13 photons s _1 cm~ 2 . The unidentified TeV source 
HESS J1813-178 (with an integrated flux above 200 GeV 
of 1.2 x 10~ 12 photons s" 1 cm" 2 ) may also be associated 
with an OB association iAharonian et al 1 l2005bfl . 

The decay of charged pions created in hadronic col- 
lisions between non-thermal and thermal ions produces a 
+ neutrino flux. This flux can be estimated from the 
7-ray flux produce d by neutral pion-decay by imp osing en- 
ergy conservation (lAlvarez-Muniz fc Halzer3l2002l) . Fig. 1231 
shows the resulting flux from models B and E, including 
neutrino oscillation. It is clear that the predicted fluxes are 
neither above the atmospheric neutrino background (ANB), 
nor above the sensitivity limit of the proposed IceCube de- 
tector. Thus, ou r finding is contrary to the conclusion of 
iBednarekl i2005l) . who predicted a detectable muon neutrino 
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Figure 23. The differential spectra of muon neutrinos produced 
by neutral pion decay from models B and E (the fluxes have been 
reduced by a factor of 2 to account for neutrino oscillation). The 
atmospheric neutrino background flux within 1° of the source 
iLiparilll993l) is s hown, toge ther with the 3-yr sensitivity of the 
IceCube detector <Hil]||200ll) . 



flux from WR 20a. This discrepancy is due to Bednarek's use 
of a higher conversion efficiency of kinetic power into non- 
thermal particles, largely through the assumption rj = 1. 



5.4 Colliding winds in dense stellar clusters 

Our work is also relevant to more complicated sys- 
tems involving colliding winds, such as dense stellar clus- 
ters. Hydrodynamical simulations of the mul tiple collid- 
ing w inds in clusters have been performed by iRaga et alJ 
q200 lTI . with particular applic ati on to the Galac ti c Cen - 
ter in ICoker fc Melial (1199^. ICoker fc Pittardl i2005t) . 
iRockefeller et alJ l)2005|h and lCuadra et alJ l)2006|h A net- 
work of bubbles forms around the stars, which are embed- 
ded in a medium of hot shocked gas. Highly non-linear DSA 
should occur at the stellar termination shocks in stellar clus- 
ters if the pre-shock B-field is smal l enough (B < 10 fiG, 
lEllison. Decourchelle fc Ballet! Eooj) . though the most ex- 
treme modification (e.g., compression ratios much greater 
than 4) will occur only at the cluster's edge. 

Shock modification is also a potential explanation for 
the lower than expected temperature of the X-ray emit- 
ting gas observed in sev eral dense stellar clusters (e.g., 
IStevens fc Hartwelll 12003) . though mass-loading and slow 
equilibration of the post-shock ion and electron tempera- 
tures are other possibilities. In addition, the X-ray lumi- 
nosities of the diffuse gas are generally higher than ex- 
pected, which indicat es mass- loading and/or a low ther- 
maliz ation efficiency dStevens fc Hartwelll 120031 : lOskinoval 
l2005lh the latter also being consistent with efficient par- 
ticle acceleration and shock modification. If the rate 
at which energy is removed from the cluster (for ex- 
ample, by radiative cooling) exceeds 30% of the en- 
ergy injection rate, a time independent cluster wind 
cann ot form llSilich. Tenorio-Tagle fc Rodriguez-Gonzalez! 
l2004h . Though the removal of energy through particle ac- 
celeration was not considered , efficiencies of this order are 
expected (e.g.. iBvkovl Il999lh in which case a full time- 
dependent solution for the cluster wind is required. 



Models of the Colliding Winds Binary WR140 25 



Amongst the dense stellar clusters known in our Galaxy, 
the central cluster around Sgr A* is unusual in many re- 
spects, including the recent detectio n of apparently th ermal 
gas with a temperature kT « 8 keV <Muno et alJl2004T) . and 
TeV emission. Observations with HESS show a 1-1 TeV 
luminosity of ~ 10 35 erg s~ 1 l|Aharonian et alJl2004h . The 
upper limit on the TeV source size is < 3' (95% confidence 
level), which corresponds to < 7 pc at the distance of the 
Galactic Center. The photon spectrum above the 165 GeV 
threshold is specified by F(E) = F E^ y , with F w 2.2 
and F « 2.5 x 1CT 12 photons s~ 1 cm~ 2 Te V" 1 , and has 
been independently confirmed with MAGIC l|Albert et alJ 
l200rJ) . Various origins of the TeV emission have been 
proposed, including the cent ral supermassive black hole 
Aharonian fc Neronavll2005al ) . the young SN R Sgr A Eas t 
Crocker et al.ll2005l) . the dark matter halo jHornsI Exjit) . 
and from non-thermal particles accelerated at t he multiple 
colliding winds shocks ijQuataert fc Loebll20f)5f ). which we 
now discuss. 

The total wind power estimated from the cluster of 
stars within ~ 10" (~ 0.4 pc) of the Galactic Center is 
~ 3 x 10 38 ergs _1 . In the context of colliding winds, the 
TeV emission may arise from IC upscattering of ambient 
photons by non-thermal electrons, or from the decay of neu- 
tral pions created by the collisions of non-thermal and ther- 
mal hadrons. We now consider each process in turn. It is 
likely that IC cooling limits the maximum Lorentz factor 
that non-thermal electrons within the ce ntral parsec can 
attain to ~ 10 jQuataert fc LoeblEooH) . which is only 
marginally high enough for IC scattering to account for the 
detected TeV flux. The TeV emission occurs from electrons 
whose cooling time is less than the flow time out of the 
cluster i.e. from those with 7 > y c . For a power- law dis- 
tribution of electrons 71(7) oc 7 _p , the IC photon spectral 
index is F = (p + 2)/2, for 7 > 7 G . Hence the TeV spec- 
tral index of 2.2 implies that p « 2.4. If instead neutral 
pions are responsible, the observed spectral index indicates 
that p ~ 2.2. In the context of multiple colliding winds, 
lOzernov. Genzel fc Usovl (jl997Tl estimates the emission from 
pion-decay to be ~ 3 x 10 36 ergs -1 , assuming 10% accelera- 
tion efficiency for hadrons at the multiple shocks, and a 10% 
energy transfer into pions i.e. a 1% efficiency overall. Other 
estimates in the literature sug gest a somewhat lower lumi- 
nosity from this process (e.g.. lQuataert fc Loebll2005f) . but 
this process nevertheless appears to be plausible on energy 
grounds. 

The required value of p for each of these processes is 
much softer than we find here for WR 140, and also for the 
theor etical limit for multiple shocks fp=1.0. fPoDe fc Melrosel 
Il994h . This implies that re-accele ration at multiple shocks 
within the cluster, as proposed bv lOzernov. Genzel fc Usovl 
l)l997fh does not occur. It is possible to imagine several rea- 
sons why this might be the case. Firstly, there is unlikely 
to be a seed population of non-thermal particles from wind- 
embedded shocks encountering the termination shock of the 
winds because of the large adiabatic cooling experienced as 
the plasma flows out to ~ 0.01 — 0.1 pc. Hence, the par- 
ticle spectrum immediately downstream of the termination 
shocks is likely to correspond to the single-shock limit. Sec- 
ondly, it is not entirely clear that there will be numerous 
opportunities for particle re-acceleration within the radius 
of the cluster, since the cooling timescale of the shocked gas 



is long, and it remains hot, and therefore subsonic. Analyti- 
cal solutions of cluster winds reveal that the shocked gas, on 
average, only becomes supersonic at the outer radius of the 
cluster iChevalier fc Clegall985fl . Within the cluster inte- 
rior, the network of hot post-shock gas experiences pressure 
gradients and weak shocks as it escapes, but in general is not 
likely to be repeatedly shocked. This lack of opportunity for 
re-acceleration means that p is likely to remain close to the 
value obtained from a single shock. 

Perhaps the most natural explanation is that the ma- 
jority of particle acceleration occurs at low Mach number 
shocks (M ~ 3 and M » 5, for IC or pion-decay emission 
respectively) within the network of hot shocked gas. This 
conclusion does not rule out the possibility of modification at 
the termination shocks of the individual winds, which would 
be consistent with the low X-ray te mperature of the c ooler 
emission component (kT ~ 0.8 keV. lMuno et al1l2004l) . but 
does exclude shock modification as the only process since 
a hard particle spectrum (p < 2) is expected at the high 
particle energies (7 > 10 J ) which are necessary. An alter- 
native explanation is sub-diffusive cross-field transport (see 
Sec. l2~Ol . 

The Galactic Center region has also been de- 
tected at GeV energies. The observed luminosity of the 
EGRET source 3EG J1745-2852 is sa 2 x lO^ergs" 1 
iMaver-Hasselwander et all Il998l : lHartman et alJ Il999f) . 
However, while its position is consistent with Sgr A*, the 
angular resolution of 1° covers a solid angle about 100 times 
larger than the emission seen with HESS. As the EGRET 
spectrum is about 4 orders of magnitude brighter than the 
observed HESS emission extrapolated as a power-law to GeV 
energies, different processes/sources may well be responsible 
for the detected emission (although different assumptions 
concerning the injection spectrum and diffusion coefficients 
of non-therm al protons may account for bo th observations, 
as shown by lAharonian fc Neronovl l2005bf) . Regardless of 
the origin of non-thermal protons, if the emission is from 
pion-decay, a substantial TeV neutrino flux should be seen 
by neutrino telescopes. 

While it remains difficult to draw firm conclusions on 
the nature of the TeV emission observed from the Galac- 
tic Center, we note that the nearby Arches cluster, which 
has an integrated mass-loss rate from its stellar components 
within an order of magnitude of the Galactic Center cluster, 
has not yet been detected at TeV energies. This may argue 
against the cluster wind scenario as an explanation for the 
TeV emission from the Galactic Center. 



6 SUMMARY AND FUTURE DIRECTIONS 

In this paper we have applied models to the observed emis- 
sion from WR 140 at <j> = 0.837, when the radio emission is 
near maximum, and when X-ray data are available. Predic- 
tions for the X-ray and 7-ray emission from IC scattering, 
relativistic bremsstrahlung, and neutral pion decay, as well 
as for the muon neutrino flux, are made. Our calculation s 
are based on the orbital solution of iDouehertv et all i2005l) , 
and contain a number of simplifying assumptions, includ- 
ing i) that the non-thermal electrons and ions are acceler- 
ated at the shocks bounding the WCR (this is assumed to 
be through DSA, but need not necessarily be so) and have 
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power-law spectra with a sharp cutoff at 7 max = fO 5 before 
advecting downstream, ii) the non-thermal particle distri- 
bution at these shocks is normalized with respect to the lo- 
cal post-shock thermal energy density, the non-thermal ions 
have 100 x the energy density of the non-thermal electrons, 
and neither the spectral index of the energy distribution 
nor 7 m ax vary along the shocks, iii) the energy density of 
the magnetic field is similarly normalized to the local post- 
shock thermal energy density, iv) particle acceleration has 
no effect on the thermal shock structure i.e. there is no shock 
"modification", v) 1C cooling occurs in the Thomson limit 
with isotropic scattering, vi) the emission from secondary 
particles is ignored. Rapid 1C cooling of the high energy 
electrons in the post-shock flow means that the local value 
of 7max is typically << 10 5 . The assumptions 7 max = 10 5 
and £ re i,i = 100£ re i,e, together with the assumptions in v), 
mean that the predicted emission at GeV energies and above 
is relatively uncertain at this stage. Our predictions for the 
VER1TAS-4 flux, in particular, should be treated as an es- 
timate. 

Despite the large number of assumptions listed above, 
it is possible to make a number of key deductions con- 
cerning the nature of WR 140 and the particle acceleration 
process at the WCR. We clearly show that VLBA imaging 
observations of WR 140 do not currently have enough 
sensitivity to constrain the wind momentum ratio 77, 
through determination of the shock opening angle. VLB1 
observations are only sensitive to high-brightness emission 
that arises in the region of the apex of the WCR, where 
the WCR has yet to attain an asympototic opening angle. 
Furthermore, we find that the relative positions of the stars 
and the peak of the observed emission from the WCR is not 
necessarily related to 77, due to occultation by the stellar 
winds. Such observations have been used to determine 77 
in a number of wider CWB systems where occulatation 
effects are expe ct ed to be relatively unimportant (e.g. 
[ Dougherty et alj Il996l: iDoug hertv. Williams fc Pollaccol 
hood IContreras et al.lll997l : IWilliams et alJll997l) . but ~ 
caution against the application of this method to closer 
systems like WR 140. In spite of these limitations in inter- 
preting VLBA observations, it is possible to constrain the 
wind momentum ratio using simultaneous fits to the radio 
and X-ray data to 77 w 0.02. A remaining concern is that the 
line-of-sight angle into the system in our preferred model 
is some way from the observationally determined value, 
although the uncertainty in the latter is not negligible. 
Additional data is required to further constrain the orbital 
parameters of WR 140. 

The stellar mass-loss rates in our preferred model are 
Mwr « 4.3 x 10~ 5 M yr _1 and Mo « 8 x 10~ 7 M Q yr _1 . 
While the O star mass-loss r ate is substant i ally l ower 
than the values estimated by iRepolusT et alj <l2004l) for 
04-5 supergiants, an order of magnitude reduction in the 
mass-loss rate is consistent with the recent estimates of 
iFullerton. Massa fc Pririial i2006t) . and with some current 
explanations of the observed near-symmetry of X-ray lines. 
Such sizeable reductions have extensive implications for 
massive star evolution (including gamma-ray bursts), and 
feedback on the interstellar medium. There is also some 
scope for reducing th e stellar luminosities adopted from 
IDoughertv 1 et all feOOEft . 

An important finding is that models of the radio data re- 



quire p < 2. Several possible mechanisms can in principle ex- 
plain this. In any case, a flatter spectrum for the non-thermal 
particles has important implications for the modelling of 
the high energy non-thermal emission from such systems, 
with almost all previous models in the literature assuming 
that p — 2. While our fits are currently non-unique, models 
where the Razin effect is responsible for the low frequency 
radio turnover (which allow p = 2) can be ruled out on the 
basis that the implied efficiency of electron acceleration is 
extraordinarily high. Furthermore, they predict the detec- 
tion of IC emission by INTEGRAL which has not yet been 
achieved, and in some cases non-thermal X-ray emission at 
fluxes higher than the observed thermal emission. Observa- 
tions with Suzaku should be more sensitive than those made 
with INTEGRAL, and will allow tighter constraints to be 
placed on the strength of the magnetic field. Our inability to 
obtain fits to the radio data with p > 2 does not necessarily 
exclude the possibility of shock modification, but on balance 
we believe that this is unlikely to occur in WR 140. How- 
ever, large non-linear effects , comparable to those possible 
at the reverse shock in SNRs llEllison. Decourchelle fc Ballet! 
120051) . may occur at the wind- wind collisions in stellar clus- 
ters where the B-field is much lower. 

Our flux predictions at orbital phase 0.837 are sum- 
marized in Tabled The photon fluxes in the INTEGRAL, 
GLAST, and VERITAS-4 band-passes range from 3x 10" 7 - 
3 x 10~ 4 , 10 -9 - 10~ 6 , and 4 x 10 -15 - 3 x 10~ 13 , respec- 
tively, with our preferred model having fluxes of 3 x 10~ 6 , 
2 x 10~ 8 , and 5 x 10~ 14 photonss -1 cm~ 2 . The high energy 
non-thermal emission from WR 140 is expected to display 
significant variability with orbital phase as the stellar sep- 
aration and line-of-sight angle change, both of which affect 
the opacity of absorption in the stellar radiation fields and 
the emission from the resulting IC cascade of secondary elec- 
trons and positrons. In addition, the spectral index and the 
acceleration efficiency of the primary non-thermal particles 
may also change. We emphasize that the future detection of 
TeV emission from CWB systems will almost certainly in- 
dicate pion-decay, since the high photon fluxes prevent the 
acceleration of electrons beyond 7 ~ 10 5 — 10 e . WR 140, 
along with WR 146 which has the brightest non-thermal ra- 
dio e mission of any CWB llDoughertv. Williams fc Pollaccol 
l200d). and WR 147 which is the widest CWB system on the 
sky (|Wilhain^e^alJ^997j) and the only system resolved in 
X-rays dPittard et al.ll2002f) . present a real opportunity for 
the unambiguous detection of pion-decay from accelerated 
ions. The estimated flux of TeV photons from CWB sys- 
tems compares favourably to some other mechanisms, and 
may contribute part or all of the flux from some TeV sources. 
The multiple colliding winds in the Galactic Center cluster 
may make a non-negligble contribution to the observed TeV 
flux. However, wide CWB systems are not expected to be 
strong neutrino sources. 

Tighter constraints on p and the nature of the shocks 
and particle acceleration in colliding winds systems (binaries 
and clusters) are best achieved through the dual strategy 
of improving the theoretical models and obtaining obser- 
vational data at MeV and GeV energies. On the theoretical 
side, models should explicitly account for the transfer of part 
of the available pre-shock energy into accelerated particles, 
self-consistently calculate the non-thermal electron and ion 
spectra and their subsequent evolution as they advect with 
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the downstream flow, and include Klein-Nishina effects and 
anistropic IC scattering. Upcoming observations by AGILE 
and GLAST at MeV and GeV energies, where the emission is 
no longer degenerate with B and p, should then distinguish 
between the non-unique model fits to the radio data, and 
determine the B-field in the WCR and the spectral index of 
the non-thermal particles. 
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